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Figure  1:  A  3-user  interference  channel. 

1  The  Studied  Problem 

Interference  is  a  major  impairment  to  successful  communication  in  commercial  and  mili¬ 
tary  wireless  systems.  In  mobile  ad  hoc  networks,  interference  is  created  when  different 
transmitters  share  the  same  channel  resource.  Uncontrolled  interference  reduces  data  rates 
throughout  the  network  and  causes  outages  at  dense  locations.  The  medium  access  control 
protocol  typically  limits  the  number  of  simultaneous  conversations  and  consequently  the 
system  performance.  Interference  is  thus  a  critical  impairment  in  ad  hoc  networks.  Com¬ 
munication  in  the  presence  of  interference  is  often  analyzed  using  an  abstraction  known 
as  the  interference  channel.  In  the  example  interference  channel  of  Fig.  1,  three  different 
transmitters  wish  to  communicate  with  three  receivers.  Each  transmitter  has  a  message 
only  for  its  paired  receiver.  Assuming  the  transmitters  share  the  same  time  and  frequency 
resources,  each  transmission  creates  interference  at  the  unintended  receivers.  There  may 
be  other  sources  of  interference,  not  illustrated,  such  as  jamming  in  military  networks,  or 
self- interference  created  from  nonlinearities  in  the  radio  frequency  components;  those  sources 
are  not  captured  in  the  basic  interference  channel. 

The  general  capacity  of  the  interference  channel  and  the  design  of  practical  schemes 
approaching  the  known  upper  bounds  on  sum  rates  have  been  of  great  interest  over  the 
last  30  years.  The  earliest  attempts  to  characterize  the  capacity  region  of  the  interference 
channel,  inspired  by  the  framework  established  by  Shannon  in  [1],  were  focused  on  two- user 
interference  channels.  Although  the  special  cases  of  strong  and  very  strong  interference  have 
been  solved  [2,3],  the  general  capacity  of  the  interference  channel  is  still  an  open  problem. 
The  difficulty  lies  in  the  dependence  between  the  receivers  such  that  the  existing  limiting 
expressions  [4,5]  require  optimization  over  many  variables  rendering  them  useless  for  practical 
purposes.  Recently,  a  series  of  attempts  have  been  made  to  describe  an  approximation  of 
the  asymptotic  sum  capacity  behavior  known  as  the  maximum  achievable  multiplexing  gain 
or  degrees  of  freedom  (DoF)  [6]  where  the  focus  is  on  the  high  SNR  regime  and  constructing 
interference- free  signals  at  the  receivers.  The  DoF  studies  paved  the  way  for  a  novel  method 
of  dealing  with  interference,  known  as  interference  alignment  [6,7]. 

Interference  alignment  is  a  linear  precoding  technique  in  which  users,  without  directly 
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Figure  2:  A  3-user  2x2  MIMO  interference  channel  where  using  IA  interference  at  each 
receiver  is  confined  to  a  single  dimension. 

talking  to  each  other,  cooperatively  and  simultaneously  transmit  signals.  In  contrast  to  other 
interference  management  techniques  such  as  orthogonal  access,  decoding  the  interference 
[4]  or  treating  the  interference  as  noise  [8],  interference  alignment  achieves  the  maximum 
multiplexing  gain  in  a  /v-user  interference  channel  [6]  such  that  nodes  sharing  the  same 
interference  channel  can  achieve  rates  equal  to  one  half  their  corresponding  interference-free 
rates. 

Consider  a  linear  system  of  equations  with  two  equations  and  three  variables.  This  sys¬ 
tem  can  represent  the  received  signal  at  two  antennas  where  the  first  variable  is  the  desired 
signal,  the  other  two  are  the  interference,  and  the  coefficients  are  the  channel  values.  In 
general,  this  underdetermined  system  does  not  have  a  unique  solution.  Assume,  however, 
that  the  coefficients  of  the  interference  messages  in  the  second  equation  differ  from  the  first 
equation  only  by  a  scale  factor.  In  that  case,  it  is  possible  to  get  a  single  equation  with 
the  desired  signal  as  the  only  variable  by  creating  a  linear  combination  of  the  two  equations 
and  solve  for  the  desired  message  interference  free.  Therefore,  for  a  general  underdeter- 
mined  system  of  equations,  if  there  exists  a  relationship  between  coefficients  of  (some  of) 
the  variables,  one  can  potentially  eliminate  a  subset  of  them  and  if  the  resulting  system  was 
proper,  a  unique  solution  for  some  of  the  original  variables  can  be  found.  This  simple  idea 
is  the  basis  of  interference  alignment  where  the  equations  are  a  combination  of  the  received 
signals  at  multiple  antennas,  times  instances,  and/or  frequency  tones  and  each  receiver  is 
only  interested  in  a  subset  of  the  variables  (the  symbols  transmitted  by  its  corresponding 
transmitter).  The  mentioned  dependency  between  the  coefficients  also  gives  a  structure  to 
the  interference. 

From  an  algebraic  point  of  view,  the  goal  of  interference  alignment  is  confining  the  ag¬ 
gregate  interference  at  each  receiver  such  that  it  occupies  a  smaller  subspace  than  it  would 
without  any  cooperation  between  the  transmitters.  Therefore,  each  transmitter,  rather  than 
optimizing  its  own  link,  cooperates  with  other  transmitters  to  manage  the  interference  it 
causes  to  non-intended  receiver.  More  technically,  interference  alignment  reduces  the  dimen¬ 
sionality  of  interference  in  the  received  signal  space  and  effectively  reduces  the  number  of 
discernible  interferers.  Fig.  2  depicts  the  same  3-user  2x2  MIMO  interference  channel  as  in 
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Fig.  1  where,  in  this  case  using  IA,  the  interference  at  each  receiver  is  confined  to  a  single 
dimension  enabling  the  system  to  transmit  3  streams  simultaneously  compared  to  orthogonal 
channel  access  where  only  2  interference-free  streams  can  be  transmitted  simultaneously. 

The  objective  of  this  project  was  to  establish  the  viability  of  interference  alignment  in  mo¬ 
bile  ad  hoc  networks  (MANETS)  under  practical  assumptions.  Several  problems  were  posed 
and  solved  that  provide  insight  into  when  and  how  interference  alignment  can  contribute  to 
overall  performance  increases  in  MANETS. 

1.1  Problems  Solved  in  Interference  Alignment 

This  project  studied  several  important  problems  in  the  domain  of  interference  alignment. 
This  section  provides  highlights  on  the  most  significant  results.  More  details  are  provided 
on  the  solution  to  each  problem  in  Section  2. 

1.1.1  User  Arrival  in  MIMO  Interference  Alignment  Networks 

Because  of  the  dynamic  nature  of  networks,  it  may  be  desirable  to  add  a  set  of  secondary 
users  who  desire  access  to  the  channel  already  being  used  by  a  set  of  user  cooperating 
through  interference  alignment.  In  our  work,  we  found  the  minimum  number  of  secondary 
transmit  antennas  required  so  that  a  secondary  user  can  use  the  channel  without  affecting 
the  sum  rate  of  the  active  users,  under  a  zero-forcing  equalization  assumption.  When  the 
secondary  users  have  enough  antennas,  we  derived  several  secondary  user  precoders  that 
approximately  maximize  the  secondary  users’  sum  rate  without  changing  the  sum  rate  of 
the  active  users.  When  the  secondary  users  do  not  have  enough  antennas,  we  performed 
numerical  optimization  to  find  secondary  user  precoders  that  cause  minimum  degradation 
to  the  sum  rate  of  the  active  users.  Through  simulations,  we  confirmed  that  i)  with  enough 
antennas  at  the  secondary  users,  gains  equivalent  to  the  case  of  all  the  users  cooperating 
through  interference  alignment  is  obtainable,  and  ii)  when  the  secondary  users  do  not  have 
enough  antennas,  large  rate  losses  at  the  active  users  can  be  avoided.  These  results  have 
appeared  in  [9]. 

1.1.2  MIMO  Interference  Alignment  in  Random  Access  Networks 

It  is  likely  that  interference  alignment  will  not  be  applied  simultaneously  across  all  users  in 
the  network.  Rather  MANETs  that  employ  interference  alignment  will  support  clusters  of 
users,  each  performance  interference  alignment.  In  our  work,  we  analyzed  a  multiple-input 
multiple-output  (MIMO)  interference  channel  where  nodes  are  randomly  distributed  on  a 
plane  as  a  spatial  Poisson  cluster  point  process.  Each  cluster  used  interference  alignment 
(IA)  to  suppress  intra-cluster  interference  but  unlike  most  work  on  IA,  we  did  not  neglect 
inter-cluster  interference.  We  also  connect  the  accuracy  of  channel  state  information  to 
the  distance  between  the  nodes,  i.e.  the  quality  of  CSI  degrades  with  increasing  distance. 
Accounting  for  the  training  and  feedback  overhead,  we  derived  the  transmission  capacity 
of  this  MIMO  IA  ad  hoc  network  and  then  compared  it  to  open-loop  (interference-blind) 
spatial  multiplexing.  Our  simulations  reveled  several  exemplary  system  setups  where  spatial 
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multiplexing  outperformed  IA  due  to  the  imperfect  channel  state  information  or  the  non- 
aligned  inter-cluster  interference.  These  results  have  been  submitted  for  publication  in  [?]; 
a  shorter  version  has  already  appeared  in  [?]. 

1.1.3  Interference  Alignment  with  Distributed  Antennas 

Distributed  antennas  are  one  way  to  provide  more  robust  transmission  in  wireless  systems. 
The  antennas  might  be  distributed  across  several  base  stations,  or  located  on  different  sides 
of  a  vehicle  for  example.  Existing  IA  solutions  cannot  be  applied  to  distributed  antenna 
systems  (DAS)  as  they  neglect  the  per-remote-radio  power  constraints  imposed  on  distributed 
precoders.  This  research  solved  the  problem  of  how  to  performance  IA  in  a  distributed 
antenna  system.  Two  different  constraints  were  considered:  ones  with  a  limit  on  maximum 
per-radio  power,  and  ones  with  a  strict  equality  constraint  on  per-radio  power.  The  rate- 
loss  incurred  by  a  simple  power  back-off  strategy,  used  in  systems  with  maximum  power 
constraints,  was  characterized  analytically.  It  was  also  shown  that  enforcing  strict  power 
constraints  avoids  such  a  rate-loss  but  negatively  affects  IA  feasibility.  For  such  systems,  an 
IA  algorithm  was  proposed  and  feasibility  conditions  are  derived  based  on  the  concept  of 
system  properness.  These  results  have  been  submitted  for  publication  in  [10]. 

1.1.4  Prototyping  MIMO  Interference  Alignment 

The  performance  of  IA  depends  on  the  practical  issues  such  as  the  performance  of  synchro¬ 
nization,  channel  estimation  and  feedback.  In  this  paper,  a  prototype  is  implemented  for  the 
IA  system  with  three  users.  There  have  been  IA  prototypes  in  recent  years,  but  the  previous 
prototypes  have  not  considered  the  distribution  of  the  nodes  in  IA  network.  The  nodes  are 
physically  distributed  in  our  prototype  not  sharing  their  time  and  frequency  references  with 
any  other,  thus  working  independently,  which  enables  the  experimental  study  of  I A  under 
the  most  practical  setup.  For  the  distributed  system,  the  over-the-air  schemes  for  time  and 
frequency  synchronization  and  analog  feedback  are  studied  and  implemented.  According  to 
the  measurement  from  our  prototype,  it  is  shown  that  IA  achieves  the  sum  rate  from  the 
previous  analysis  on  imperfect  channel  information.  In  addition,  some  other  measurements 
are  performed  considering  the  accuracy  of  IA  solution,  synchronization  of  nodes,  and  CSI 
feedback.  For  the  accuracy  of  IA  solution,  the  performance  of  IA  versus  the  number  of  itera¬ 
tions  for  an  iterative  IA  method  is  measured.  For  synchronization  accuracy,  the  performance 
with  different  residual  frequency  offset  is  measured.  Finally,  for  the  CSI  feedback  quality, 
analog  feedback  and  scaler  quantization-based  limited  feedback  is  first  compared.  These 
results  are  about  to  be  submitted  for  publication  in  [?].  Earlier  related  results  were  reported 
in  [11], 

2  Summary  of  the  Key  Results 

2.1  User  Arrival  in  MIMO  Interference  Alignment  Networks 

In  general,  a  high  computational  complexity  [12, 13]  and  overhead  [14, 15]  cost  is  associated 
with  finding  the  IA  precoders  and  combiners.  For  example,  overhead  is  required  to  acquire 
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propagation  channel  state  information,  exchange  channel  state  information,  and  to  coor¬ 
dinate  the  transmissions  and  reception  of  multiple  users.  This  overhead  is  incurred  often, 
essentially  every  time  the  channel  changes.  Further,  the  amount  of  overhead  scales  super  lin¬ 
early  with  the  number  of  users.  The  premise  of  this  project  is  applied  locally  in  small  closed 
clusters  to  compromise  between  increased  data  rate  and  moderate  overhead.  This  premise 
was  established  in  prior  work  by  the  PI,  where  he  established  the  optimality  of  performing 
interference  alignment  over  small  groups  of  users  in  a  variety  of  settings  including  overhead 
and  with  channel  estimation  error. 

Once  a  set  of  users  have  aligned  their  interference,  it  is  desirable  to  retain  their  alignment 
status  until  the  channels  change.  An  example  of  a  set  of  existing  IA  users  is  illustrated  in 
Fig.  3.  Most  future  interference-limited  wireless  networks,  however,  will  be  packet-switched 
with  bursty  data  traffic,  requiring  frequent  changes  in  the  number  of  active  users  [16].  Also, 
the  feedback  overhead  [17]  and  the  number  of  antennas  at  each  node  [18]  practically  limit 
the  number  of  user  pairs  that  can  cooperate  through  IA.  Therefore,  in  an  IA  network,  there 
will  be  nodes  that  are  arriving  to  the  network  and,  although  not  included  in  the  existing  IA 
setup,  wish  to  communicate  with  their  receivers.  An  example  of  new  nodes  joining  an  IA 
network  is  depicted  in  Fig.  3.  Determining  when  such  nodes  can  be  admitted  to  the  network 
(are  allowed  to  transmit)  and/or  developing  transmission  strategies  for  them  is  referred  to  as 
user  arrival.  Without  developing  appropriate  user  arrival  techniques  for  the  new  users,  IA 
can  not  be  utilized  in  a  dynamic  network  without  incurring  overhead  to  reconfigure  the  entire 
network.  This  motivates  considering  the  problem  of  user  arrival  in  interference  alignment 
networks. 

In  this  work,  we  consider  a  set  of  active  users  exhausting  the  network  resources  by 
cooperatively  utilizing  IA  and  a  further  set  of  secondary  users  who  wish  to  communicate  in 
this  network.  Such  a  network  topology  can  be  a  small  part  of  a  larger  network  which  is,  due 
to  the  limited  number  of  antennas  at  each  node  and  physical  distance,  clustered  into  separate 
regions.  Fig.  4  depicts  an  ad  hoc  network  clustered  into  separate  IA  groups  and  some  nodes 
which  are  not  part  of  any  I A  cluster.  We  assume  the  secondary  users  are  required  to  have 
minimum  impact  on  the  performance  of  the  active  users,  defined  as  the  sum  rate  of  the  active 
users,  and  the  active  users  ignore  the  presence  of  the  secondary  users  when  designing  their 
precoders  and/or  equalizers. 

The  main  technical  contributions  of  this  work  are  summarized  as  follows.  We  compute 
a  zero-impact  threshold  for  the  number  of  secondary  transmitter  antennas  where  secondary 
transmitters  with  more  (or  equal)  antennas  than  this  threshold  can  use  the  communications 
medium  without  degrading  the  sum  rate  of  the  active  users.  We  find  optimum  and  sub¬ 
optimum  secondary  user  precoders  for  two  cases:  (i)  when  the  secondary  users  satisfy  the 
zero-impact  threshold  and  (ii)  when  they  do  not. 

When  the  secondary  users  satisfy  the  zero-impact  threshold,  there  exists  a  set  of  precoders 
that  will  not  degrade  the  sum  rate  of  the  active  users.  Thus,  the  secondary  users  can  optimize 
an  objective  function  of  their  own  fink  by  selecting  a  precoder  from  this  set.  We  choose  the 
achievable  sum  rate  as  the  objective  function  and  derive  optimum  and  suboptimum  precoders 
for  the  case  of  one  and  two  secondary  users  respectively.  For  more  than  two  secondary 
users,  we  choose  the  degrees  of  freedom  (DOF)  as  the  objective  function,  defined  as  the 
slope  of  average  sum  rate  (b/s/Hz)  versus  logarithm  of  signal-to-noise  ratio  (dB)  at  high 
transmit  power.  We  propose  successive  IA  precoding  and  show  it  is  optimum  for  various 
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Figure  3:  A  set  of  existing  IA  users  fully  utilizing  network  resources  and  a  set  of  new  nodes 
joining  the  network. 

network  setups  determined  by  the  number  of  active  users,  secondary  users,  and  antennas  at 
each  node.  When  the  secondary  users  have  fewer  antennas  than  the  zero-impact  threshold, 
we  search  for  secondary  user  precoders  causing  minimum  degradation  to  the  sum  rate  of 
the  active  users  through  a  steepest  descent  search  over  the  Grassmann  manifold.  For  this 
numerical  optimization,  we  propose  three  initial  solutions  of  varying  degrees  of  complexity. 

2.1.1  Mathematical  Framework  and  Problem  Statement 

Consider  a  iFa-user  MIMO  interference  channel  where  the  ?'th  transmitter  and  receiver  are 
equipped  with  M%  and  Nt  antennas,  respectively.  In  general,  each  transmitter  i  uses  a 
precoding  matrix  F,  of  dimension  Mi  x  di  to  transmit  di  streams  to  its  corresponding  receiver; 
which  is  then  decoded  by  the  ith  receiver  after  processing  the  received  signal  with  a  combining 
matrix  W*.  As  precoders  with  orthonormal  columns  are  preferred  in  MIMO  channels  [19], 
we  assume  Fj  £  O  where  O  is  the  set  of  matrices  with  orthonormal  columns.  For  each  time 
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Figure  4:  An  ad  hoc  network  partitioned  into  separate  I A  clusters  which  some  trans¬ 
mit/receive  pairs  are  not  part  of  any  I A  clusters. 


instant,  the  received  signal  at  the  zth  receiver,  with  perfect  timing  and  synchronization,  is 


Ka 


y*  =  ^2  HifcFfcXfc  +  Zj  i  =  l,...,K 


k= 1 


(i) 


where  H^.  is  the  matrix  of  channel  coefficients  of  a  block  fading  channel  between  transmitter  k 
and  receiver  i,  the  transmitted  signal  from  the  Ah  node  is  x,  with  power  constraint  E{x*x(}  = 
P  and  z i  is  the  AWGN  with  elements  in  CJ\f  (0.  a2)  where  a2  is  the  noise  power  spectral 
density  which  includes  thermal  noise  and  white  excess  interference  from  unaccounted  sources. 
Moreover,  we  also  assume  H ^  has  full  column  rank  for  all  i,k  implying  enough  scattering 
in  the  channel. 

The  goal  of  IA  for  the  MIMO  constant  channel  is  to  design  the  precoding  and  combining 
matrices  such  that  the  following  conditions  are  satisfied 


f  rank(WiHiiFi)  =  d; 

\  WjHjfcFfc  =  0  Vk^i 


(2) 


We  assume  the  system  of  IA  in  (2)  is  feasible  [18].  Feasibility  of  IA  implies  that  inter¬ 
ference  from  the  active  transmitters  at  the  Ah  active  receiver  is  confined  to  an  Nt  —  (fa 
dimensional  subspace;  assume  columns  of  C(  are  a  non-unique  basis  for  this  interference 
subspace.  We  further  assume  each  active  receiver  uses  a  zero-forcing  equalizer,  given  by 
Wf  =  [Idi,0]  [HjjFj,  Cj]_1.  It  can  be  shown  that  the  total  achievable  sum-rate  of  the 
active  users  is  given  by 


R 


au 

sum 


Ka  d, 


z^Z^log 

2=1  n=  1 


To/  di 


((F-H. 


„PiH«F 


n 


(3) 


where  P,  =  (Dq  —  C*C*)  is  the  projection  matrix  into  the  nullspace  of  the  interference 
subspace  at  the  Ah  active  receiver,  q0  =  ^  and  en  is  the  nth  row  of  In . 

Assume  Ks  secondary  users  request  sharing  the  network  resources.  Define  K  =  Ka  +  Ks, 
and  without  loss  of  generality,  let  secondary  users  be  the  last  Ks  users  in  the  ordered  set  of 


user  indices  /C  =  {1, . . . ,  Xa,  +  1, . . . ,  X}.  We  assume  the  first  Ka  users  do  not  change 
their  precoders  or  receiving  filters  when  the  secondary  users  join  the  network.  This  has 
two  important  implications:  a)  the  active  transmitters  do  not  help  the  transmission  of  the 
secondary  users,  b )  the  active  receivers  are  not  aware  of  the  interference  from  the  secondary 
users. 

In  the  presence  of  the  secondary  users,  R^m  changes  from  (3)  to 


Ka  di 

EE 

i—1  n=  1 


log 


lo/di 


F*  TET  *  T3  U  I? 


)“1+Ef= 


ica+i  |WIH,FfcF*H*W* 


(4) 


where  we  have  assumed  that  H,^.  is  independent  of  Hm„  for  i  ^  n  or  k  ^  rn.  We  divide 
the  contribution  of  this  work  into  three  parts  corresponding  to  answering  the  following  three 
questions 


1.  How  many  antennas  does  a  secondary  node  need,  Ms,  in  order  to  use  the  network  re¬ 
sources  without  degrading  the  active  user’s  sumrate  ((4)  is  equal  to  (3))?  Equivalently 


Ms\Mi  >  Ms 


i  >  Ka  — >■ 


K 

E  =  0 

-TS  II  ^ 


(5) 


2.  When  the  secondary  nodes  have  enough  antennas,  what  is  their  optimum  precoding 
matrices  such  that  in  addition  to  not  degrading  the  sumrate  of  the  active  users  the 
achievable  sumrate  of  the  secondary  nodes  is  also  maximized?  Equivalently,  we  solve 

k  /  K 

argmax  V  logdet  ( I  +  +  V 

rw<-r  «eatJZ+i  V  d‘  ,-u*kd' 

(6) 

s.t.  [lPiHI1;)' . iPa,HAi1)TFi.-0  k=Ka  +  l,...,K. 


3.  When  the  secondary  nodes  do  not  have  enough  antennas,  what  is  their  optimum  pre¬ 
coding  matrices  such  that  the  sum-rate  degradation  at  the  active  users  is  minimized. 
Equivalently,  we  solve 


argmax  ££  log  f  1  +  )  .  (7) 

r'‘»« . V  e»  (Qi  +  gW,  (Eh&+i  H^FtFJHi)  W*j  e*J 

where  Q,  =  (F'H'PiHaF,)-1. 


2.1.2  Threshold  on  the  Number  of  Secondary  Transmitter’s  Antennas 

We  now  present  a  summary  of  the  key  results  obtained  at  this  work.  We  first  present  the 
threshold  on  the  minimum  number  of  transmit  antennas  required  at  the  secondary  nodes 
such  that  secondary  nodes  can  communicate  without  degrading  the  sum-rate  of  active  nodes. 
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Then  we  present  the  results  on  two  scenarios,  one  when  the  secondary  transmitters  have 
enough  antennas  (more  than  this  threshold)  and  the  other  when  the  secondary  transmitters 
are  lacking  enough  number  of  transmit  antennas.  With  enough  antennas,  the  secondary 
transmitters  can  avoid  causing  any  sum-rate  degradation  at  the  active  receivers  where  we 
derive  optimum  secondary  precoders  that,  in  addition,  maximize  the  achievable  sum-rate  of 
the  secondary  nodes.  Without  enough  antennas,  the  secondary  transmitters  are  bound  to 
decrease  the  sum-rate  of  the  active  nodes  and  in  this  case  we  present  secondary  transmitter 
precoder  matrices  that  minimize  this  sum-rate  loss. 

As  discussed  before,  each  active  receiver  uses  a  zero-forcing  to  receiver  to  combine  its  sig¬ 
nal.  The  zero-forcing  receiver  is  based  on  projecting  the  received  signal  into  the  nullspace  of 
the  interference  subspace  at  each  receiver.  Therefore,  if  the  interference  from  the  secondary 
transmitters  is  also  aligned  into  this  interference  subspace,  it  will  be  cancelled  without  chang¬ 
ing  the  signal-to-noise-plus-interference  (SINR)  at  the  active  receivers.  Using  this  fact,  in  [20] 
we  prove  the  following  Lemma. 

Lemma  1  The  kth  >Ka  transmitter  can  send  dk  streams  without  degrading  (3)  if 

Ka 

Mk  >  ^  ]  di  +  dk  k  £  {Ka  +  1, . . . ,  K}.  (8) 

1=1 

Lemma  1  states  that  if  no  sum-rate  degradation  at  the  active  nodes  is  tolerable,  the  number 
of  secondary  transmit  antennas  need  to  be  at  least  equal  to  the  sum  of  the  streams  in  the 
active  network,  du  phis  number  of  the  desired  streams  at  the  secondary  pair  itself. 

Note  that  the  only  parameter  from  the  active  network  is  the  total  number  of  streams  in  that 
network  and  the  obtained  threshold  is  independent  of  the  number  of  antennas  at  the  active 
nodes. 

2.1.3  Enough  Antennas  at  the  Secondary  Transmitters 

The  exact  solution  of  (6)  is  out  of  the  scope  of  this  work  (it  is  equivalent  to  solving  the 
sum-rate  maximizing  solution  of  the  MIMO  interference  channel)  and  we  focus  on  several 
special  case. 

2.1.4  Single  Secondary  Pair 

Let  the  columns  of  V*,  G  O  span  the  right  nullspace  of  H/  =  [(PxHife)’  V  •  •  >  (P Ka^-Kak)*]*  ■ 
In  [20]  we  prove  that  for  a  single  secondary  pair,  a  solution  for  (6),  Fxa+ i,  is  given  by  Lemma 
2. 

Lemma  2  For  Ks  =  1,  FT  sovling  (6)  equals  V^Ga,  where  the  columns  of  Gk  are  the  ds 

K 

most  significant  eigenvectors  of  V]*H]T,(Ijvfc  +  ^  and  k  =  Ka  +  l. 

i  l .//  a-  1 


Two  Secondary  Pairs  When  Ka  1  and  Ks  =  2,  in  [20]  we  approximately  solve  (6) 
using  Corollary  1.  Similar  to  Lemma  2,  let  Fa,  =  V^Ga,.  Then  we  have 
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Corollary  1  When  Ka  1,  (6)  is  approximately  solved  by  setting  the  columns  of  G&  to  the 
ds  most  significant  eigenvectors  of 


Basically,  in  [20]  we  show  that  when  Ka  1,  the  interference  of  the  other  secondary  trans¬ 
mitter  at  each  secondary  receiver  becomes  negligible  and  the  problem  folds  back  to  the  case 
of  Ks  =  1. 


More  Than  Two  Secondary  Pairs  Solving  (6)  for  Ks  >  2  is  equivalent  to  solving 
the  general  sum-rate  maximizing  precoder  of  the  MIMO  interference  channel  and,  to  date,  a 
closed-form  solution  directly  solving  (6)  for  Ks  >  2  does  not  exist  [21].  We  provide,  however, 
a  precoder  design  called  successive  interference  alignment  [20]  which  maximizes  the  multi¬ 
plexing  gain  of  the  secondary  network  (also  called  the  achievable  degrees-of-freedom)  for  a 
certain  (and  plausible)  network  configuration.  In  successive  interference  alignment,  the  sec¬ 
ondary  transmitters  first  align  their  interference  at  the  active  receivers  into  the  corresponding 
interference  subspaces.  Then,  the  secondary  transmitters  and  the  secondary  receiver  coop¬ 
eratively  modify  the  secondary  precoding  matrices  such  that  a  new  interference  alignment 
network  is  formed  within  the  secondary  nodes.  We  can  only,  however,  provide  a  conjecture 
on  the  optimality  of  successive  interference  alignment  in  certain  network  configurations  and 
its  performance  evaluation  is  left  for  future  work. 

Conjecture  1  Consider  a  3K-user  interference  channel  for  K  G  Z+  where  the  transmit¬ 
ter/receiver  pairs  are  divided  into  K  groups  of  3  users  each,  {G\,  G2 , . . . ,  Gk},  such  that  the 
nodes  of  the  ith  group  have  3i  —  1  antennas.  Performing  successive  IA  on  the  kth  group, 
2  <  k  <  K ,  through  creating  effective  channels  between  the  nodes  of  Gk  based  on  the  inter¬ 
ference  subspaces  and  precoders  of  the  {1, . . . ,  k  —  1}  groups  achieves  the  same  DoF  as  if  all 
the  3 K  users  had  done  IA  together. 


2.1.5  Not  Enough  Antennas  at  the  Secondary  Transmitters 

In  general,  (7)  is  a  complex  non-convex  problem  and  its  exact  solution  is  out  of  the  scope 
of  this  contribution.  By  assuming  Ks  =  1,  we  simplify  (7)  so  that  a  numerical  optimization 
technique  can  be  utilized  to  search  for  the  optimum  secondary  precoder.  The  simplified 
objective  problem,  however,  is  still  non-convex  and  the  better  the  initial  solution  provided 
to  the  numerical  optimization  problem  the  better  the  final  solution  of  the  numerical  search. 
Next  we  present  three  initial  solutions  with  varying  degrees  of  complexity. 

Alternating  Minimization  Based  Initial  Solution  Let  W  =  0^  0^=1  v^enWjHjs, 
s  =  0[=i  0n=i  ^enQ,;e*  ,  and  d  =  Kada.  It  is  inferred  from  (2)  that  S  is  full  rank  and  at 
high  SINR,  (7)  can  be  written  as 

argmin  logdet  (i  +  -^(Ij®  F*)W*S_1W(Id~  0  Fs)  )  .  (9) 

Fsee>  \  da  J 


11 


Let  the  eigenvalue  decomposition  of  W*S-1W  be  If  the  rows  of  Ij  0  F*  were 

equal  to  a  linear  combination  of  eigenvectors  corresponding  to  zero  eigenvalues  in  the 
det  in  (9)  would  attain  its  minimum  value  of  1.  Therefore,  ideally 

I<J  0  Fs  =  LJ-yy-As,  (10) 

where  the  columns  of  are  the  columns  of  corresponding  to  zero  eigenvalues  of 
W*S-1W  and  As  G  O  is  an  (Ms  —  1)  di  x  (Ylf= i  di)ds  combining  matrix. 

In  [20],  we  show  that  (10)  does  not  have  an  exact  solution  and  we  present  an  iterative 
algorithm  to  approximately  solve  (10)  in  the  least  squares  sense.  Assume  Fs  in  (10)  is  given. 
We  seek  an  As  minimizing  |  |I^  <S>  Fs  —  U-vpAs|||.  Equivalently 

A-g  =  argmax  91c  (tr  (A*U*y  (Id~  0  Fs)^  ,  (11) 

where  9Te{ . }  select  the  real  part  of  a  complex  value.  We  solve  the  optimization  problem  of 
(11)  using  the  solution  to  the  “Procrustes  problem”  [22],  Now  consider  (10)  but  assume  that 
this  time  As  is  given.  We  seek 

Fs  =  argmin  |  |Ij0  Fs  -  U^yAs |||  =  argmax  91c  (tr  ((Ij  0  FS)A*U^))  .  (12) 

Fsee>  Fsee>  v  v  J  J 

Similar  to  (11),  in  [20]  we  solve  (12)  using  the  general  solution  to  the  Procrustes  problem. 

Interference  Leakage  Minimizing  Initial  Solution  Compared  to  the  alternating  mini¬ 
mization  solution,  a  less  complex  (and  less  accurate)  initial  solution  can  be  found  by  revisiting 
the  constraint  on  the  secondary  user  precoders  required  when  no  sum  rate  degradation  at 
the  active  receivers  is  tolerable.  Instead  of  requiring  the  secondary  transmitters’  interfer¬ 
ence  to  be  confined  in  the  interference  subspace  of  each  active  receiver,  we  minimize  (in 
the  least  squares  sense)  the  interference  leakage  caused  by  the  secondary  users  by  solving 
argmiiip  eC,  ||HSFS|||.  A  solution  to  this  problem  is  given  by  setting  the  columns  of  Fs  to 
the  ds  least  significant  right  singular  vectors  of  Hs. 

DoF-Preserving  Initial  Solution  In  both  the  previously  presented  initial  solutions  in¬ 
terference  from  the  secondary  transmitter  is  not  confined  to  the  interference  subspaces  of 
the  active  receivers.  Thus,  the  degrees-of-freedom  of  the  active  users’  network  is  zero.  By 
aligning  the  secondary  transmitter’s  interference  at  some  of  the  active  receivers,  however,  a 
non-zero  degrees-of-freedom  can  be  achieved.  Thus,  we  seek  a  solution  to  the  problem  of 
argminFs eC)  |HSF.,|  |0.  Instead  of  directly  solving  this  combinatorial  problem  [23],  we  find 
the  largest  subset  of  active  receivers  such  that  interference  from  the  secondary  transmitter 
can  be  perfectly  aligned  at  the  active  receivers  and  by  constructing  F,s  based  on  this  subset 
of  active  receivers  a  non-zero  multiplexing  gain  is  achieved  at  the  active  network.  Note  that 
this  solution  is  optimal  for  da  =  1. 
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Figure  5:  Traditional  single  cluster  analysis  (left)  versus  analyzing  a  large  network  with 
many  such  clusters  in  it  (right). 


2.2  MIMO  Interference  Alignment  in  Random  Access  Networks 

In  large  networks,  such  as  MANETs,  I A  will  be  used  independently  in  separate  clusters 
and  the  nearby  nodes  that  are  not  coordinating  with  any  one  cluster  will  cause  non-aligned 
interference  at  the  receivers  (an  example  of  this  configuration  is  shown  Fig.  4).  This  is 
motivated  by  the  following  observations. 

1.  The  number  of  antennas  at  each  node  is  a  limiting  factor.  It  is  shown  in  [24]  that  the 
number  of  nodes  that  can  cooperate  through  MIMO  IA  is  limited  by  the  number  of 
antennas  at  each  node. 

2.  Overhead  practically  limits  the  cluster  size.  The  overhead  of  IA  grows  super  linearly 
with  the  number  of  users  [25],  and  hence  it  is  likely  that  small  groups  of  nodes  will 
coordinate  to  perform  IA. 

3.  More  cooperation  is  not  always  better.  Recently,  [26]  showed  that  because  of  inherent 
channel  uncertainty,  there  is  a  moderate  cluster  size  above  which  spectral  efficiency  at 
best  saturates,  and  in  many  practical  scenarios  (e.g.  when  pilots  are  used  for  channel 
estimation),  actually  decreases  if  more  nodes  join  the  cluster  to  cooperate. 

In  the  context  of  aggregate  large  network  performance,  single  cluster  analysis  (e.g.  DoF 
studies)  does  not  capture  the  impact  of  interference  from  the  other  nodes  in  the  network 
and  can  lead  to  unrealistic  cooperation  gains  which  would  not  be  attainable  if  the  inter- 
cluster  interference  was  accounted  for  [26,27].  An  example  of  single  cluster  analyzes  versus 
analyzing  a  large  network  with  many  such  clusters  is  shown  in  Fig.  5 

Therefore,  the  premise  of  this  project  is  that  IA  is  applied  locally  in  small  closed  groups 
dividing  a  large  network  into  separate  clusters.  This  premise  was  also  established  in  prior 
work  by  the  PI,  where  he  established  the  optimality  of  performing  interference  alignment  over 
small  groups  of  users  in  a  variety  of  settings  including  overhead  and  with  channel  estimation 
error.  When  dealing  with  large  networks,  a  relevant  metric  of  the  system  performance  is 
the  transmission  capacity  [28],  defined  as  the  number  of  successful  transmission  per  unit 
area,  subject  to  a  constraint  on  outage  probability.  The  transmission  capacity,  in  contrast 
to  other  network-wide  system  performance  metrics  such  as  transport  capacity,  generally 
leads  to  closed-form  expressions  or  tight  bounds  providing  insight  into  the  network  design 
parameters  [29] .  Also,  transmission  capacity  analysis  generally  take  into  account  large-scale 
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Figure  6:  An  Aloha- like  cluster- wise  channel  access  strategy. 

fading  effects,  such  as  pathloss,  which  is  also  a  key  factor  in  determining  the  benefits  of 
interference  alignment  in  practical  scenarios  as  for  example  cooperation  between  some  nodes 
that  are  not  in  each  other’s  transmission  range  or  are  not  causing  significant  interference 
to  their  respective  receivers  may  not  be  necessary  and  even  can  be  counterproductive.  This 
motivates  finding  the  transmission  capacity  of  decentralized  networks  employing  cluster-wise 
interference  alignment. 

In  this  work,  we  consider  a  large  ad  hoc  network  where  nodes  are  partitioned  into  separate 
clusters  each  cooperating  through  IA.  We  assume  a  four-stage  transmission  protocol.  In  the 
first  stage,  with  a  finite  length  training  period ,  imperfect  CSI  for  the  cross  links  is  obtained 
through  MMSE  channel  estimation.  In  the  second  stage,  the  estimated  CSI  is  fed  back 
to  the  other  nodes  in  the  cluster  during  the  feedback  period.  In  the  third  stage,  the  IA 
transmit/receive  filters  are  computed.  In  the  last  stage,  using  the  rest  of  the  finite- length 
channel  block,  the  nodes  communicate  data  using  a  cluster-wise  slotted  Aloha-like  channel 
access  protocol  where  at  random,  all  nodes  in  a  cluster  either  transmit  simultaneously  or  turn 
off  their  transmission.  Such  a  topology  can  represent  a  MANET  where  cooperation  between 
close-by  nodes  is  obtainable  and  the  disjoint  clusters  can  synchronize  their  transmission  using 
GPS  or  low-overhead  message  passing.  An  example  of  Aloha-like  cluster-wise  channel  access 
is  shown  in  Fig.  6. 

The  main  technical  contributions  of  this  work  are  summarized  as  follows.  We  first  derive 
the  point-to-point  probability  of  outage  in  such  a  network  accounting  for  the  imperfect 
channel  estimation.  Then,  we  derive  the  optimum  training  period  assuming  a  quasi-static 
channel  maximizing  cluster  throughput  while  taking  into  account  the  feedback  overhead. 
Finally,  we  derive  the  transmission  capacity  of  this  clustered  ad  hoc  network  and  compare 
it  with  a  much  simpler  parallel  system  where  a  single  node  in  each  cluster  uses  spatial 
multiplexing.  The  goal  is  gaining  insight  into  the  operation  regimes  where  IA  outperforms 
common  transmission  techniques  by  considering  node  density,  channel  block  length,  transmit 
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Figure  7:  An  instance  of  the  transmitter’s  distribution  when  the  number  of  transmitters  per 
cluster  is  K  =  3. 

power,  number  of  available  antennas,  and  the  number  of  nodes  in  each  cluster. 

2.2.1  Mathematical  Framework  and  Problem  Statement 

The  spatial  locations  of  the  potential  transmitters,  $,  are  modeled  as  a  planar  Neyman- 
Scott  cluster  point  process  [30].  In  this  process,  the  cluster  centers  are  modeled  by  a  parent 
homogeneous  Poisson  point  process  (PPP)  <E>p  of  density  Ap.  Each  parent  point  x  G  <Fp 
forms  the  center  of  a  cluster  around  which  K  daughter  points  are  uniformly  distributed  in 
a  circle  of  radius  R.  The  resulting  process1  is  a  stationary  point  process  of  density  K Ap. 
We  assume  clusters  randomly  access  the  channel  with  probability  Pa  effectively  reducing  the 
density  of  this  PPP  to  Ap  =  PaAp.  The  receiver  of  a  transmitter  at  x  is  denoted  by  x  and 
is  assumed  to  be  randomly  located  at  distance  DT  from  its  transmitter  forming  an  iV  x  ]V 
MIMO  link.  The  receivers  are  not  part  of  the  point  process  <E>.  An  instance  of  the  nodes’ 
location  is  shown  in  Fig.  7. 

In  this  work,  a  typical  transmitter  (a  transmitter  chosen  at  random)  is  considered  and 
its  performance  is  analyzed.  This  transmitter  is  typical  in  the  sense  that  the  performance  of 
I A  in  this  node  is  a  representative  of  the  average  performance  of  I A  in  the  network  [30,31]. 
Since  the  underlying  point  process  is  stationary,  without  loss  of  generality,  it  can  be  assumed 
that  the  typical  transmitter  is  at  the  origin.  Denote  the  cluster  to  which  the  transmitter  at 

1The  parent  points  <f>p  will  not  be  a  part  of  the  final  point  process. 
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the  origin  belongs  by  Tv  The  received  signal  at  receiver  x,  x  £  d/0,  is 


y x  ^  \J 9 x z H x 2 F z s~  +  Xc  +  u£,  (13) 

ze^o 

where  Xc  =  ]C2c<y/^,,  yhAH.;2F,sr,  is  the  inter-cluster  interference,  g%z  and  H,~  represent 
the  pathloss  and  the  matrix  of  channel  coefficients  between  the  transmitter  z  and  the  receiver 
x,  Fz  is  the  precoder  at  transmitter  z  with  the  transmitted  signal  s-,  such  that  E{s*s2}  =  P, 
and  u,  ~  CJ\f(0,  N0T)  is  the  additive  white  Gaussian  noise.  In  this  work,  it  is  assumed 
that  F*FZ  =  I,  because  of  tractability  and  the  observation  that  the  gain  attained  otherwise, 
such  as  with  the  MMSE  algorithm  in  [32]  or  the  Max-SINR  algorithm  in  [12],  is  limited 
and  confined  to  the  low  SNR  regime  where  the  inter-cluster  interference  is  generally  not 
dominant.  In  every  cluster,  channel  state  information  is  estimated  at  the  receivers  as  in  [33] 
and  conveyed  to  all  other  nodes  of  the  cluster  using  an  error-free  instantaneous  feedback 
link.  We  propose  to  model  the  uncertainty  in  the  MIMO  channels  using  a  Gauss-Markov 
model  of  the  form  [34, 35] 


H„  = 


2  TTW 

xznxz 


fixZ^XZ  %•)  Z  G 


(14) 


where  is  the  estimated  channel,  Exz  represents  the  estimation  error  with  i.i.d.  terms 
distributed  as  CAf(0, 1),  and  B2Z  is  the  normalized  variance  of  the  estimation  error.  It  is 
assumed  that  the  channel  is  quasi-static  block-fading  such  that  H  is  constant  for  a  block 
duration  of  length  T  and  then  changes  independently.  Training,  feedback,  and  data  trans¬ 
mission  are  assumed  to  be  all  orthogonal  in  time,  in  the  same  coherence  time  or  frame  T  [36]. 
Hence,  B±z  is  set  to  be  related  to  the  average  received  SNR  at  each  link,  qy, ,  as  [33,  Section 
II. B] 


B2  = 

r  XZ 


1  +  7/ lx 


2  _|_  yfoQxz 


N 


(15) 


where  70  =  X-  and  Tt  >  KN  is  the  number  of  channel  instances  spent  for  training 
H £z  [37].  For  analytical  tractability,  it  is  also  assumed  that  H"  is  used  to  construct  the 
precoders/equalizers  and  nodes  effectively  ignore  the  imperfection  in  CSI  in  their  design. 

At  each  cluster,  a  A'- user  system  of  IA  is  feasible  if  there  exists  a  set  of  matrices  W  = 
{Wz\ z  e  \i>o}  such  that,  given  the  received  signal  of  (13),  the  following  constraints  are 
met  [6]: 


f  rank  (W4H4xFx)  =  Ns 

\  W£H£zFz  =  0 


(16) 


where  is  the  combining  filter  used  at  receiver  x  and  Ns  is  the  number  of  interference- free 
streams  each  transmitter  can  send  to  its  receiver.  The  linear  equalizer  presented  in  [38] 
is  an  examples  of  a  possible  receive  filter  in  (16).  It  is  assumed  that  the  IA  precoders 
are  designed  using  the  alternating  minimization  algorithm  in  [32,  Section  III. A]  such  that 
Fx  is  independent  of  for  all  £  G  \l/.  Also,  it  is  assumed  that  the  set  of  {N,  Ns,  K} 
constitutes  a  feasible  IA  system,  which  for  the  MIMO  interference  channel  requires  that 
2N  —  (K  +  1  )Na  >  0. 
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From  (16),  interference  at  receiver  x  is  confined  to  an  TV  —  Ns  dimensional  subspace. 
Let  [{•}]  represent  horizontal  concatenation  of  the  elements  in  {•}.  Then,  as  I A  pre- 
coders/equalizers  are  constructed  using  H™  as  given  by  (14),  the  N  x  (K  —  1  )NS  matrix 
Of  J*  =  [{H”F,M  tf0}]  spans  an  N  —  Ns  dimensional  subspace.  Let  the  singular 

value  decomposition  of  J*  be  UjsSjiVj.  and  let  the  rows  of  Wx  be  the  columns  of  Ujr 
corresponding  to  zero  singular  values  in  Sj4 .  As  Wx  is  independent  of  H"’/;F,,,,  it  satisfies  the 
conditions  in  (16)  and  is  a  valid  zero-forcing  (ZF)  equalizer  for  IA.  Using  this  ZF  receiver, 
the  post-processing  SINR  of  the  nth  stream  at  receiver  x  is 
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xz^-xz 


(17) 


ze®/&0 

\ _ 


where  for  all  zGf,  =  (e* W^H^F^)*,  e^z  =  (e* W^E^F^)*,  and  en  is  the  nth  column  of 
an  Ns  x  Ns  identity  matrix.  Note  that  Is  represents  the  residual  error  from  the  direct  fink 
and,  as  the  distance  between  the  transmitter  and  receiver  is  constant,  its  pathloss  (and  the 
error  variance)  are  not  random  variables  and  so  it  is  separated  from  Ie  to  emphasize  this 
point.  Let  the  entries  of  H,~  and  E„:2  be  i.i.d.  Gaussian  terms.  As  W;,:  and  Fz  are  unitary 
matrices  independent  of  H  c  -  and  E:rz ,  due  to  the  doubly  unitarily  invariance  of  the  Gaussian 
distribution,  h,  ,  and  exz  will  be  column  vectors  of  length  Ns  with  i.i.d.  Gaussian  terms  and 
independent  of  each  other  (as  H %z  and  Exz  are  independent  of  each  other).  Note  that  (17) 
is  in  fact  independent  of  the  stream  index  n.  We  divide  the  contribution  of  this  work  into 
three  parts  corresponding  to  answering  the  following  three  questions 

1.  Given  the  SINR  expression  in  (17),  what  is  the  point-to-point  probability  of  outage? 
Equivalently,  given  the  SINR  threshold  of  9,  we  seek  to  find 

P*AW=P!°(7^>0), 

where  b  is  the  receiver  corresponding  to  a  typical  transmitter  located  at  the  origin. 
If  point-to-point  outage  is  of  interest,  this  analysis  reveals  the  dependency  of  perfor¬ 
mance  on  important  system  parameters  such  as  training  period  length,  transmit  power, 
number  of  antennas,  density  of  the  nodes,  and  pathloss. 

2.  For  a  given  block  fading  of  length  T,  what  is  the  optimum  training  length  Tt  locally 
maximizing  the  goodput  of  each  cluster?  Equivalently,  we  solve 

-  T  —  K2N  —  T 

Tt  =  argmax  - — - -KNsPlA(6)  log2(l  +  9) 

Tt  2 

s.t.  KN  <  Tt  <  T  -  K2N 

Finding  the  optimum  training  period  gives  insight  into  the  overhead  associated  with 
I A  compared  to  the  other  well  known  transmission  techniques. 
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3.  Let  q(Xp)  =  1  —  PgA($)  =  e  be  the  maximum  cluster  density  corresponding  to  the 
outage  threshold  of  9.  Accounting  for  overhead,  we  seek  the  normalized  transmission 
capacity  defined  is 


C(e) 


t-k2n 

T~ 


Tt 


q-\e)K(l 


e). 


The  transmission  capacity  hints  at  maximum  spectral  efficiency  attainable  by  MIMO 
IA  in  a  decentralized  network. 


2.2.2  Point-to-Point  Outage  Probability 

Now  we  present  a  summary  of  the  key  results  obtained  at  this  work.  In  this  section  we 
derive  the  exact  point-to-point  outage  probability  at  a  typical  receiver.  Then,  assuming 
fixed  feedback  overhead,  we  solve  for  the  optimum  training  period  locally  maximizing  each 
cluster’s  goodput.  Next,  we  derive  the  exact  transmission  capacity  for  the  case  of  a  single 
data  stream  from  each  transmitter.  Finally,  we  present  an  upper  bound  on  the  transmission 
capacity  for  the  case  of  arbitrary  number  of  streams  from  each  transmitter.  Detailed  results 
are  presented  in  [39-41,  Chapter  3]. 

In  (17),  since  h*zhcz  and  are  i.i.d.  T(NS,  1)  random  variables,  we  denote  them 

both  by  hxz  for  notational  simplicity.  Denote  the  typical  transmitter  at  the  origin  by  o. 
Considering  a  transmitter  at  the  origin  is  equivalent  to  conditioning  on  the  existence  of  a 
point  at  the  origin.  Since  every  point  belongs  to  some  cluster,  conditioning  on  the  existence 
of  a  point  at  the  origin  equals  the  presence  of  a  cluster  with  a  daughter  point  at  the  origin. 
Since  the  parent  point  process  is  a  PPP,  an  additional  cluster  with  a  daughter  point  at  the 
origin  can  be  added  to  it  without  changing  the  statistics  of  the  other  points  of  the  process. 
Succinctly,  the  Palm  probability  of  a  Newman-Scott  cluster  process  is  P°  =  P  *  \l/0  where 
*  denotes  superposition  [30].  This  implies  that  assuming  a  point  of  the  cluster  process  at 
the  origin  equals  the  original  point  process  <E>  plus  an  additional  cluster  which  has  a  point 
at  the  origin.  Also  this  additional  cluster  at  the  origin  \F0  is  independent  of  the  original 
process  *F.  Since  the  tagged  transmitter  at  origin  does  not  contribute  to  the  interference 
at  the  receiver,  it  is  convenient  to  use  the  reduced  Palm  probability  denoted  by  P!o  instead 
of  Palm  probability.  Reduced  Palm  probability  is  similar  to  Palm  probability,  except  that 
the  point  at  the  origin  is  not  considered  in  the  computation  of  the  probability  and  hence 
P!o  =  P*  {\F0\{o}}.  From  (17),  the  probability  of  success  is  therefore  PgA(6>)  =  P!o(7,5A  >  0), 
where  9  is  the  SINR  threshold. 


Theorem  1  For  the  system  model  described  in  Section  2.2.1,  the  success  probability  when 
each  cluster  uses  IA  is 
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where  Oj  {s)  and  /^(s)  are  Laplace  transforms  of  intra  and  inter-cluster  interference  given 
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where  rj  = 
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(20) 


2.2.3  Optimum  Training  Period 

For  a  given  block  fading  of  length  T,  the  transmitters  spend  Tt  >  KN  channel  instances  for 
training  the  links.  We  also  assume  a  prefect  analog  feedback  link  where  the  receivers  send 
the  trained  channels  over  a  period  of  Tf  =  K2N  channel  instances  to  the  transmitters.  In 
practice,  the  transmitters  select  Tt  to  optimize  some  performance  criteria.  In  this  work,  we 
assume  transmitters  use  the  goodput  at  each  cluster  defined  as 

T  —  K2N  —  T 

T,  =  argmax  - - - log2(l  +  6) 

Tt  J-  (21) 

s.t.  KN  <Tt  <T  —  K2N, 


where  acc0unts  for  the  transmission  opportunities  lost  due  to  overhead,  KNS  is 

the  total  number  of  streams  in  each  cluster,  and  PgA(d)  log2(l  +  6)  is  the  rate  multiplied  by 
the  times  the  connection  exists,  i.e.  SINR  passes  the  threshold  9.  Note  that  in  (21),  PgA(d) 
is  implicitly  a  function  of  the  training  period  Tt. 

For  a  given  node  mobility  and  hence  a  Doppler  frequency  f,j  ~  P,  the  training  period 
can  be  written  as  a  fraction  of  the  total  block  length  T ,  i.e.  Tt  =  5T  =  SR.  Therefore,  the 
optimization  problem  of  (21)  can  be  rewritten  as 


5  =  argmax  (1  -  5  -  fdK2N)KNsPlA{9)  log2(l  +  9) 

5 

s.t.  fjKN  <  <S  <  (1  -  UK2N), 


(22) 


where 
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d y  ,  and  £Ii(s)  is  given 


by  (20). 

With  a  convex  relaxation  on  Tt  (and  therefore  5)  to  change  its  domain  to  the  real  num¬ 
bers,  the  optimization  problem  of  (22)  is  convex  and  solvable  with  any  of  the  numerical  op¬ 
timization  algorithms  [42],  Note  that,  although  complicated,  the  derivative  of  the  objective 
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function  in  (22)  w.r.t  h  is  computable  and  evaluating  the  objective  function  or  its  derivatives 
for  any  set  of  values  is  possible.  Next  we  provide  approximate  closed  form  solutions  for  some 
common  cases. 


When  Ns  =  1,  P*A  in  (22)  simplifies  to 


p;A(V  =  1)  =  e 


+  fdDCT 


t?  +  fd(<i  +  A") 


fd)C  iM- 


(23) 


The  highly  non-linear  dependency  of  (23)  on  h  can  be  converted  into  a  polynomial  one 
following  the  Taylor  expansion  approximation  method  proposed  in  [43].  Let  g(8,  fa)  be  the 
objective  function  in  (22).  Rewrite  g(8,  ff)  as  its  Taylor  series  around  frj  =  0  (infinite  block 
length  and  hence  perfect  training)  keeping  all  the  other  variables  constant 


^  c)  Q  1  Q 

5  ~  arg  max  g (h,  fd  =  0)  +  —  |  fd=0fd  +  -  — 2  I  fd=o  fl 
5  Old  *  OJd 

s.t.  fdKN  <  6  <  (1  -  fdK2N). 


(24) 


Simplifying  the  terms  in  (24)  and  removing  the  constant  scaling  coefficients  from  the  opti¬ 
mization  problem  yields 


hi  «  argmax  —  8  +  C\  -  +  C2 

s  8  8Z  (25) 

s.t.  fdKN  <8  <  (l-fdK2N), 

where  C\  and  Co  are  given  in  Appendix  [40,  Appendix  D],  Let  hi  be  the  relevant  root  of  the 
first  derivative  of  (25).  The  optimum  training  period  will  be  given  by 

ftt  1  =  min  (max (KN,  [hxT])  ,T  —  K2N)  .  (26) 

2.2.4  Transmission  Capacity 

Let  q(Xp)  =  1  —  PgA(6l)  =  e.  Accounting  for  overhead,  the  normalized  transmission  capacity 
is 


T  —  K 2  TV  —  T 

C(e)  = - - - V^Ml  -  c).  (27) 

With  the  probability  of  successful  transmission  as  given  by  (18),  the  exact  expression  of 
g-1(e)  for  general  Ns  is  not  analytically  tractable.  Next  we  give  its  exact  expression  for  the 
case  of  Ns  =  1  and  provide  a  bound  for  general  Ns  >  1. 


Single  stream  from  each  transmitter  Using  (18) 


q{ Ap)  =  1  -  e-^ 


Zl7o  1  net 
N  ' 


Ttfo 

N 


4' :Sfi)  £ilv)  =  u 


(28) 
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where  fj  =  f)\T  f  .  Substituting  (19)  and  (20)  into  (28)  gives 


loge 


e  '  t-q  -jv 


1 +D ?  1 


1_€  ^+v+d?*r2 
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Q  (e)  =  K  = 


B(o,R ) 


ttR2  /  ( 

B(o,RV  1  N 


Hvr+lk-y-oll 


v+^r+Wx-y-dW 


Ns 

da; 


-I  K-ls 


d  y. 


fl¬ 


irt?  f 

B(o,R) 


11^— y||c 


f)+\\x-y\\' 


Ns 


dx 


K 


d  y 


(29) 


where  Tt  is  the  optimum  training  period  obtained  earlier.  Using  (29),  (27)  can  be  computed. 


Greater  than  one  stream  from  each  transmitter  When  f)  >  (Ns  —  l)/e,  using  [40, 
Lemma  2] 


(-lN)  >  1  - 


Ns- 1  „A. 

1]_p-(fi-k/e)^ZJffp-{v-k/e)goo^0hoo  r'o 


k= 0 


k\ 


A°e  (V  ~  k/e)£u  (V  ~  k/e) .  (30) 


Now,  if  fj  k/e  for  all  k  E  {0, . . . ,  Ns  —  1},  fj  —  k/e  ~  fj  and  hence  (30)  equals 
l(\)  >  1  -  (  £  e-0^^Cl{rj)CIt(rj) 

\  k= 0  K'  ) 

„  .  1  —  e 

=>  CiM  > 


^  \  — 


Eto1  i)  e-^-E 
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dx 
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(31) 


2.2.5  Intuition  Obtained  from  the  Results 

Our  results  indicate  that  selecting  IA  as  the  transmission  technique  of  choice  in  a  MANET 
should  be  based  on  the  node  density,  the  mobility  of  the  nodes,  the  transmit  power,  and  the 
characteristics  of  the  underlying  communications  medium.  For  example,  in  dense  networks 
with  high  transmit  power,  spatial  multiplexing  (SM)  over  an  orthogonal  channel  access 
strategy  such  as  time  division  multiple  access  (TDMA)  can  outperform  IA  due  to  lower 
inter-cluster  interference.  Also,  the  signal-to-noise-ratio  (SNR)  switching  point  between  IA 
and  TDMA+SM  decreases  with  increasing  density  and  mobility.  In  conclusion,  MIMO  IA 
is  best  suited  for  low  Doppler  disperse  clusters  of  nodes  capable  of  local  coordination. 

2.3  Interference  Alignment  with  Distributed  Antennas 

Interference  alignment  is  a  transmission  strategy  that  is  designed  to  maximize  the  number 
of  non-interfering  symbols  that  can  be  simultaneously  communicated  over  an  interference 
network  [6].  By  doing  so,  IA  achieves  the  maximum  degrees  of  freedom  in  a  variety  of 
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single  or  multi-antenna  interference  channels  and  consequently  allows  systems  to  approach 
capacity  in  the  high  signal-to-noise  ratio  (SNR)  regime.  In  the  more  practically  relevant 
regime  of  low-to-medium  SNR,  however,  IA’s  sum-rate  performance  has  been  shown  to  be 
suboptimal  [44-48] .  This  limits  the  applicability  of  IA  in  practical  network  with  reasonable 
transmit  power  budgets  since  it  does  little  to  improve  the  rates  for  the  disadvantaged  low- 
SNR  users  [49].  To  overcome  this  practical  shortcoming,  prior  work  has  focused  on  the 
development  of  enhanced  precoding  strategies  that  improve  on  IA’s  low-SNR  performance. 
The  unifying  concept  behind  the  algorithms  in  [47,48,50,51]  is  to  potentially  forgo  perfect 
alignment  and  optimize  objective  functions  that  are  more  tightly  related  to  system  sum- 
rate.  In  addition  to  better  algorithms,  however,  IA’s  utility  can  be  improved  by  considering 
network  architectures  that  may  be  more  suitable  for  such  a  high-SNR  transmission  strategy. 
One  such  network  is  systems  with  distributed  antennas.  In  fact,  the  merits  of  combining 
distributed  antennas  with  cooperative,  or  multi-user  transmission  strategies,  have  previously 
been  established  by  the  PI  in  the  case  of  the  multi-user  MIMO  broadcast  channel  [52,53]. 

In  this  work  we  consider  the  application  of  interference  alignment  to  networks  with  dis¬ 
tributed  antenna  transmitters.  The  motivation  for  combining  IA  with  DAS  is  two-fold.  First, 
DAS  may  help  overcome  or  avoid  IA’s  low-SNR  weakness  by  boosting  signal  power  for  distant 
users.  Second,  I A  may  constitute  an  effective  co-channel  interference  management  strategy 
for  DAS.  Since  existing  IA  solutions  neglect  per-RRU  (or  per-antenna)  power  constraints, 
and  are  thus  not  applicable  to  DAS,  we  focus  on  reevaluating  the  possibility  and  perfor¬ 
mance  of  IA  in  these  more  tightly  constrained  systems.  We  note  that  while  the  antenna’s 
geographic  separation  in  DAS  necessitates  the  consideration  of  such  power  constraints,  per- 
antenna  constraints  are  in  fact  of  practical  relevance  in  all  wireless  transceivers  that  must 
operate  with  power-efficiency  in  mind,  i.e.,  even  those  with  co-located  antennas  [54-58]. 

We  consider  I A  in  two  types  of  DAS,  ones  in  which  individual  RRU  powers  are  upper 
bounded  (called  maximum  power  constraints )  and  others  in  which  all  RRUs  must  transmit 
at  a  constant  power  level  (called  strict  power  constraints )  [56].  This  work  thus  completes  the 
preliminary  results  presented  in  [59]  which  focused  on  the  special  case  of  IA  with  per-antenna 
power  constraints  in  co-located  antenna  systems.  To  satisfy  maximum  power  constraints,  we 
consider  a  simple  strategy  of  transmit  power  back-off  and  show  that  the  back-off  procedure 
incurs  a  systematic  loss  in  sum-rate.  To  gain  a  better  quantitative  understanding  of  the 
sum-rate  loss  due  to  back-off,  we  give  an  analytical  characterization  of  the  back-off  factor’s 
statistics  under  the  simplifying  assumption  of  channels  with  equal  pathloss,  i.e.,  in  the  case 
of  co-located  antennas.  The  development  of  more  sophisticated  IA  strategies  that  satisfy 
maximum  power  constraints  and  avoid  the  loss  incurred  by  back-off  is  an  interesting  topic 
for  future  work.  In  systems  where  RRUs  must  transmit  at  a  fixed  power  level,  we  show  that 
the  addition  of  such  strict  power  constraints  negatively  affects  the  feasibility  of  IA.  In  other 
words,  realizing  IA  with  strict  power  constraints  requires  more  antennas  at  the  transmitters 
or  receivers.  To  examine  this  reduction  in  feasibility,  we  leverage  the  methodology  in  [60]  to 
derive  properness  conditions  for  I A  systems  with  strict  per-RRU  power  constraints.  While 
properness  and  feasibility  are  in  general  not  equivalent,  the  results  of  [61-63]  indicate  that 
proper  systems  are  most  often  feasible  except  in  a  few  corner  cases.  As  a  result,  properness 
can  provide  a  simple  and  sufficiently  accurate  predictor  of  IA  feasibility.  To  further  demon¬ 
strate  the  true  feasibility  of  IA  in  proper  systems,  we  present  a  simple  iterative  IA  algorithm 
based  on  the  alternating  minimization  solution  in  [51,64],  While  more  sophisticated  algo- 
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Figure  8:  K- User  MIMO  interference  channel  model  with  distributed  antenna  transmitters. 
Each  transmitter  consists  of  Arru  remote  radio  units  each  with  A^/AGrtj  antennas. 

rithms  are  an  interesting  topic  for  future  work,  the  proposed  solution  is  shown  to  work  well 
in  simulation.  Namely,  the  proposed  algorithm  avoids  the  sum-rate  loss  of  power  back-off 
and  interestingly  achieves  the  same  performance  as  unconstrained  IA. 

2.3.1  Mathematical  Framework  and  Problem  Statement 

Consider  the  AT-user  interference  channel  shown  in  Fig.  8  in  which  transmitter  k  communi¬ 
cates  with  its  targeted  receiver  k  and  interferes  with  all  other  receivers  £  ^  k.  We  assume 
that  the  system  is  symmetric,  meaning  that  the  number  of  transmit  antennas  Nt,  receive 
antennas  Nr,  and  data  streams  Ns  is  the  same  for  all  transmitter-receiver  pairs.  We  denote 
such  symmetric  systems  as  (Nt  x  Nr.  NS)K .  In  our  DAS  setup,  we  assume  that  each  trans¬ 
mitter’s  antennas  are  distributed  among  AAr.u  remote  radios,  each  with  Nt / AAnu  antennas 
as  shown  in  Fig.  8. 

Let  us  define  x*,  to  be  the  A^  x  1  symbol  vector  transmitted  by  node  k  with  E  [xfcx’ji]  =  Tvs, 
F/,.  to  be  the  Nt  x  Ns  precoding  matrix  used  by  transmitter  k,  Hy  to  be  the  Nr  x  Nt  channel 
from  transmitter  £  to  receiver  k,  and  z*,  to  be  Nr  x  1  vector  of  i.i.d  complex  Gaussian  noise 
observed  at  receiver  k  with  covariance  matrix  <J2Ijvr.  Assuming  perfect  time  and  frequency 
synchronization,  user  Ar’s  received  signal,  y  & ■  can  be  written  as 

Yfc  =  ^  +  zk.  (32) 

e^k 

To  enable  the  calculation  of  the  I A  precoders,  however,  the  channels  are  assumed  to  be 
known  perfectly  to  both  transmitters  and  receivers.  In  practice,  channel  knowledge  can  be 
obtained  by  reciprocity  or  feedback  [65,66]. 

The  precoders  Fp.  are  constrained  to  satisfy  a  total  power  constraint  given  by  |  F 1 1  p  =  P, 
where  P  is  the  total  transmit  power  available  at  each  transmitter.  Total  available  power, 
however,  is  not  the  only  constraint  on  the  transmit  precoders  F .  Practical  power  amplifiers 
and  transceiver  architectures  typically  place  a  limit  on  the  power  radiated  by  individual 
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antennas  or  by  groups  of  antennas  located  at  the  same  RRU.  To  handle  per-RRU  power 
constraints  mathematically,  we  partition  the  precoders  F/0  into  ATru;  transmit  subfilters  Fk> 

^  =  [Fp,  f!2)*,  . . . ,  Ff RRu)*l  *  as  illustrated  in  Fig.  8. 


each  of  size 


Nt 


x  Ns  such  that 


Constraints  are  then  placed  on  the  Frobenius  norm  of  each  subfilter  F[r) .  In  cases  where  per- 
antenna  power  is  restricted,  constraints  are  simply  placed  on  the  norm  of  individual  rows 
of  the  precoder  F*,  (equivalently  RRUs  can  be  thought  of  as  having  a  single  antenna 
each).  With  this  notation,  we  define  the  two  power  constraints: 


1.  Maximum  Power  Constraints  where  the  maximum  power  radiated  by  an  RRU  or  an 

antenna  is  constrained,  i.e.,  |F[r)  ||A  <  in  the  case  of  per-RRU  constraints  and 

||f^r^  || 2  <  in  the  case  of  per-antenna  constraints. 

2.  Strict  Power  Constraints  where  a  strict  equality  constraint  is  placed  on  the  power 
radiated  by  different  RRUs  or  antennas,  i.e.,  ||F^  |||.  =  in  the  case  of  per-RRU 

constraints  or  ||f^r')|||  =  jr  in  the  case  of  per-antenna  constraints. 

Sections  2.3.2  and  2.3.3  analyze  the  effects  of  these  more  stringent  constraints  on  I A  feasibility 
and  performance.  When  IA’s  performance  is  considered,  the  main  metric  of  interest  is  the 
average  sum-rate  achieved  with  complex  Gaussian  signaling  and  treating  interference  as 
noise.  Under  these  assumptions,  sum-rate  is  given  by  [67] 


K 

/  \  -1 

X]l0S2 

1,/Vr  + 

+  y  H„F»F,*HJ(  HttFuFJH^ 

k= 1 

V  l^k  ) 

(33) 


where  denotes  expectation  over  the  distributions  of  the  channels  H k£  V/.\  L 

While  IA  can  be  used  with  any  receiver  architecture,  alignment  can  be  intuitively  under¬ 
stood  by  examining  the  operation  of  a  linear  interference  zero-forcing  receiver.  Define  W{; 
to  be  the  Nr  x  Ns  linear  zero- forcing  combiner  used  by  receiver  k.  The  received  signal  at  the 
output  of  Wfc  is  given  by 

W*kyk  =  W*HfcfcFfcxfc  +  W*  HmF,x*  +  W*kzk.  (34) 

i+k 

At  the  output  of  Wfc,  the  conditions  that  ensure  perfect  alignment  can  be  stated  as 


W^H^F^  =  0Ns, 
rank  (W^H^F*,)  =  Ns, 


Vfc,£e/C,  (35) 

Vfc  e  JC,  (36) 


where  /C  is  the  set  of  users  and  Cfys  is  the  Ns  x  Ns  all-zeros  matrix.  Interference  alignment 
and  cancellation  is  guaranteed  by  (35),  and  (36)  ensures  the  decodability  of  user  k' s  desired 
signal.  The  IA  conditions  in  (35)- (36)  have  been  used  extensively  in  the  literature  to  derive 
several  I A  algorithms  [50,51,64]  and  various  performance  results  [60,62,65,66,68]. 

The  work  in  [60]  leveraged  conditions  (35)  and  (36)  to  characterize  the  systems  in  which 
IA  is  feasible.  The  authors  of  [60]  defined  the  notion  of  proper  IA  systems  in  which  the 
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number  of  free  variables  in  the  transmit  precoders  F/c  always  exceeded  the  number  of  con¬ 
straints  imposed  by  condition  (35).  It  was  then  shown  that  system  properness  constitutes 
a  necessary  condition  for  I A  feasibility.  More  recent  work  has  shown  that,  while  properness 
does  not  rigorously  guarantee  IA  feasibility,  proper  systems  are  most  often  feasible  except 
in  a  few  corner  cases  [61,62],  As  a  result,  the  notion  of  properness  can  provide  a  sufficiently 
accurate  predictor  of  IA  feasibility.  For  the  symmetric  systems  considered  in  this  paper,  for 
example,  an  IA  system  is  considered  proper  (and  most  likely  feasible)  as  long  as 

Nt  +  NT>(K  +  l)Na.  (37) 

The  conditions  in  (35)- (36)  have  similarly  been  used  to  derive  iterative  algorithms  in  cases 
where  closed  form  solutions  do  not  exist.  In  the  alternating  minimization  solution  of  [51], 
the  total  power  of  leakage  interference,  defined  as  ||WfcH^F^||2,  is  iteratively  minimized 
over  alternating  choices  of  F^  and  W*,.  Using  the  derived  algorithms,  I A  was  ultimately 
shown  to  provide  good  high-SNR  sum-rate  in  a  variety  of  MIMO  interference  channels. 

Since  existing  I A  algorithms,  feasibility  results,  and  performance  analysis  do  not  consider 
per-RRU  or  per-antenna  power  constraints,  it  remains  unclear  if  IA’s  promise  carries  over 
to  DAS.  To  see  this,  note  that  algorithms  that  neglect  per-RRU  power  constraints  yield 
precoders  with  an  unbalanced  power  profile  across  different  antennas.  This  power  imbalance 
can  be  significant  in  DAS  where  different  antennas  experience  significantly  different  pathloss. 
If  maximum  power  is  constrained,  transmitters  can  back-off  their  total  transmit  power  to  en¬ 
sure  that  no  RRU  exceeds  its  power  constraint.  Power  back-off,  however,  will  result  in  a  loss 
of  both  effective  SNR  and  IA  sum-rate  which  we  characterize  analytically  in  Section  2.3.2.  In 
systems  that  require  all  RRUs  to  transmit  at  the  same  power,  existing  IA  algorithms  cannot 
be  used  altogether  and  power  back-off  cannot  be  used  to  balance  RRU  transmit  power.  In 
fact,  Section  2.3.3  shows  that  placing  such  strict  constraints  on  per-RRU  power  significantly 
affects  IA  feasibility.  For  such  systems,  Section  2.3.3  derives  revised  properness  conditions 
and  provides  an  iterative  algorithm  that  enables  calculating  constrained  I A  precoders. 

2.3.2  IA  Sum  Rate  loss  with  Maximum  Power  Constraints 

When  a  limit  is  placed  on  per-RRU  transmit  power,  IA  precoders  must  satisfy  the  following 
AAru  inequality  constraints 

||Fir)|||<-^,  Vre77,  Vfce/C,  (38) 

tVRRU 

where  7 Z  is  the  set  of  RRUs.  Since  per-antenna  constraints  are  mathematically  equivalent  to 
(38)  with  Arru  =  Ay,  we  focus  on  the  general  case  of  per-RRU  constraints.  Examining  the 
IA  conditions  given  by  (35)  and  (36),  we  note  that  the  constraints  in  (38)  should  not  affect 
IA  feasibility.  Indeed,  suppose  that  there  exists  a  set  of  precoders  F*,  and  combiners  W/,: 
that  satisfy  (35)-(36)  with  a  maximum  RRU  transmit  power  of  [3k  =  ma yireiz  HFjU’H^.  In 

this  case,  the  scaled  precoders  v  y>^v'l'  l;"  ~pk  simultaneously  ensure  that  alignment  is  preserved 
and  transmit  power  constraints  are  respected. 

While  the  constraints  in  (38)  do  not  affect  feasibility,  satisfying  them  via  power  back-off 

causes  a  systematic  reduction  in  IA  sum-rate  since  v/~P^'RRU  <  1  with  probability  one.  To 
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gain  a  quantitative  understanding  of  the  performance  degradation  induced  by  power  back¬ 
off,  we  examine  the  mean  loss  in  sum-rate  at  high  SNR  which  we  define  in  the  following 
proposition. 


Lemma  3 

is  given  by 


The  mean  loss  in  sum-rate  resulting  from  transmit  power  back-off  at  high  SNR 
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Rioss  —  NSE 


k= 1 


log2  f 


NrrvPI 


(39) 


Given  Lemma  3,  all  that  remains  to  complete  the  characterization  of  Rioss  is  to  derive 
the  statistics  of  the  random  variables  The  distribution  of  jff.,  however,  is  tied  to  the 
distribution  of  the  precoders  F/;  and  is  thus  dependent  (non-trivially)  on  the  statistical 
model  used  for  the  propagation  channels  H^.  This  makes  deriving  the  distribution  of  fff 
for  general  channel  models  intractable.  Thus,  to  simplify  the  rate  analysis,  we  make  the 
following  channel  assumption. 


Assumption  1  We  assume  that  all  channels  are  Rayleigh  fading,  i.e.,  have  i.i.dCJ\f(  0, 1) 
entries. 


Mathematically,  Assumption  1  enables  us  to  use  the  following  result  from  [68]. 

Theorem  2  Assuming  Rayleigh- fading  channels,  the  precoders  F*,  G  CNtxNs  for  k  G  K 
generated  using  the  IA  algorithms  in  [6,  51,  64]  are  Haar- distributed,  i.e.,  they  are  uniformly 
distributed  over  the  set  of  orthogonal  Ns- frames  in  CNt . 

Theorem  2  facilitates  the  derivation  of  f3f  s  statistics  by  (i)  identifying  a  single  tractable 
distribution  for  the  precoders  F*,,  and  (ii)  consequently  indicating  that  all  variables  /A  are 
statistically  equivalent.  As  a  result,  we  henceforth  drop  the  subscript  (•)*,  from  f3k- 

For  brevity,  we  give  the  final  results  on  the  distribution  of  the  power  back-off  factor  (3 
and  refer  the  reader  to  [69]  for  more  detailed  derivations. 

Proposition  1  In  the  case  of  single  stream  transmission,  i.e.,  Ns  =  1,  and  assuming  that 
is  a  constant  natural  number,  the  CDF  of  the  random  variable  f32  can  be  approximated 
as 

/  N  \Nt 

{(32  <  x}  «  <2  Ntx-  -±-  .  (40) 

V  -‘Vjuuj  / 

and  the  approximation  can  be  proven  to  be  exact  in  the  limit  of  Nt  — >  oo. 

Proposition  2  The  distribution  of  the  maximum  RRU  transmit  power  (3 2  in  the  case  of 
Ns  >  1  can  be  approximated  as 

P/32  {f2  <  x}  ~  Q  (NtX;  *•  (41) 

V  -‘Vjuuj  / 

and  the  approximation  can  again  be  proven  to  be  exact  in  the  limit  of  Nt  — >  oo. 

Having  found  the  distribution  of  the  power  back-off  factor  (3 2  for  all  Ns  >  1,  evaluating  the 
mean  loss  in  sum-rate  give  in  Proposition  3  reduces  to  a  simple  one-dimensional  integral 
that  can  be  evaluated  numerically.  Numerical  examples  about  the  loss  expected  from  power 
back  off  can  be  found  in  [69]. 
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2.3.3  IA  Feasibility  with  Strict  Power  Constraints 

To  ensure  that  systems  transmit  at  full  power,  strict  equality  constraints  can  be  placed  on 
the  IA  precoders  F&.  Recalling  the  definition  in  Section  2.3.1,  such  constraints  can  be  written 
as 

IIFf’llt  =  trace  (F<r)F<r>*)  =  A-,  Vr  €  R.Vfc  €  K.  (42) 

v  '  tVRRU 

As  stated  in  Section  2.3.2,  if  per-RRU  transmit  power  is  neglected,  IA  algorithms  will  gen¬ 
erate  precoders  in  which  the  quantities  ||F^||i?  are  continuous  random  variables.  Thus  for 
existing  IA  solutions  ||F^||^  ^  nUw  with  probability  one.  Unlike  maximum  power  con¬ 
straints,  however,  (42)  cannot  be  satisfied  via  simple  transmit  power  back-off.  Further,  it 
is  possible  that  enforcing  the  strict  power  constraints  in  (42)  will  fundamentally  affect  the 
feasibility  of  IA  in  DAS  and  will  necessitate  the  development  of  improved  IA  algorithms  that 
explicitly  account  for  per-RRU  power. 

To  explore  the  reduction  in  feasibility  caused  by  enforcing  strict  power  constraints,  we 
leverage  and  extend  the  notion  of  properness  which  was  introduced  in  [60] .  Since  the  MIMO 
interference  channels  considered  in  this  paper  are  symmetric,  properness  can  be  determined 
by  examining  the  number  of  free  variables  and  unsatisfied  equations  in  the  full  set  of  IA 
conditions  (35),  i.e.,  there  is  no  need  to  consider  subsets  of  equations  from  (35)  [60].  Let  Nv 
be  the  total  number  of  free  variables  in  the  IA  system,  and  let  and  n!(2>  be  the  number 
of  non-trivially  satisfied  equations  in  (35)  and  (42)  respectively. 

From  [60],  the  number  of  free  variables  Nv  and  unsatisfied  equations  in  (35)  are 
given  by 


Nv  =  K(Nt  +  Nv-2Ns)Ns,  (43) 

lVe(1)  =  K[K  -  1  )NS2.  (44) 

The  details  of  this  counting  argument  are  provided  in  [60].  In  short,  Nv  is  determined  by 
counting  the  free  variables  in  Ffc  VA:  e  /C  and  W/,.  VA:  e  /C  after  reducing  them  to  their  unique 
basis  representations,  while  Ne1'1  is  determined  by  counting  the  scalar  equations  in  (35).  To 
count  the  number  of  non-trivially  satisfied  equations  in  (42),  we  provide  the  following  result. 

Proposition  3  The  number  of  non-trivially  satisfied  equations  in  (42)  is  given  by 

Aj2)  =  max{A(iVRRU  -  As),0}  .  (45) 


Combining  (43),  (44)  and  (45),  the  following  characterization  of  system  properness  for  IA 
with  strict  per-RRU  power  constraints  can  be  obtained. 

Theorem  3  A  symmetric  I  A  system  with  strict  per-RRU  power  constraints  is  proper,  and 
is  thus  expected  to  be  feasible,  if  and  only  if 

(Nt  +  Nt)Ns  >(K  +  1  )NS2  +  max  {(1VRRU  -  Ns),  0}  (46) 
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K=  3,  N=2,  Single  Antenna  RRUs 
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K=  4,  N  =3,  Two-Antenna  RRUs 
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Figure  9:  Tables  highlighting  network  configurations  in  which  IA  is  infeasible,  strictly  feasi¬ 
ble,  or  feasible  only  without  per-RRU  or  per-antenna  power  constraints  as  predicted  by  the 
derived  conditions  on  system  properness. 


Examining  the  properness  conditions  derived  in  Proposition  3  for  IA  with  strict  per-RRU 
power  constraints,  we  make  the  following  observations: 

1.  For  single  antenna  RRU’s,  the  properness  condition  in  (46),  simplifies  to  the  per- 
antenna  constrained  case  examined  in  [59].  Namely,  since  Arru  =  At  >  As,  the 
condition  for  properness  simplifies  to 

(Ar  +  Nt)Ns  >  (K  +  1  )AS2  +  (At  -  Ns).  (47) 

2.  Interestingly,  per-RRU  power  constraints  only  reduce  feasibility  in  cases  where  Ay.ru;  > 
As.  In  cases  where  Ay.Ru  <  Ay.  the  properness  condition  in  (46)  reduces  to  the 
traditional  IA  properness  condition  in  (37).  As  a  special  case,  this  regime  includes 
traditional  co-located  antenna  systems  (Arru  —  1)  with  only  a  total  transmit  power 
constraint. 

3.  Per-RRU  power  constraints  destroy  the  symmetry  of  alignment  since  the  properness 
condition  no  longer  depends  on  At  and  Ay  only  through  their  sum  Ay  +  Ay.  Therefore, 
unlike  in  the  traditional  IA  case  explored  in  [60],  if  an  (At  x  Ar,  NS)K  system  is  proper 
and  thus  likely  feasible,  the  ((At  +  1)  x  (Ar  —  1),NS)K  system  need  not  be  feasible. 
The  reciprocal  (Ar  x  At,  NS)K  system  need  not  be  feasible  either. 

4.  Interestingly,  in  the  case  of  As  =  1  with  single  antenna  RRUs,  i.e.,  Ay.Ru  =  At,  transmit 
antennas  are  entirely  useless  for  alignment  or  any  transmit-side  interference  nulling. 
This  can  be  seen  by  noting  that,  in  this  case,  the  properness  condition  simplifies  to 
(At  +  Ar)  >  (A  +  1)  +  (At  —  1)  =t-  Ar  >  K.  The  condition  Ar  >  K  implies 
that  the  receivers  must  have  enough  antennas  to  single-handedly  cancel  all  unaligned 
interference  and  decode  their  desired  signals.  Note  that  in  this  case,  the  transmit 
precoders  are  simply  equal-gain  beamforming  vectors  [54]  which  we  have  now  shown 
can  never  be  used  for  alignment. 
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5.  The  condition  in  (46)  confirms  the  intuition  that  having  multiple  antenna  RRUs  (such 
that  iVRR u  <  At)  significantly  reduces  the  effect  of  per-RRU  power  constraints  and 
thus  improves  IA  feasibility. 

To  get  a  better  understanding  of  the  effect  of  per-RRU  power  constraints  on  IA  feasibility, 
Fig.  9  tabulates  some  example  scenarios  in  which  IA  is  always  feasible,  feasible  only  in  the 
absence  of  per-RRU  power  constraints,  or  always  infeasible,.  Comparing  the  two  K  =  4  cases 
with  single  and  multi-antenna  RRUs,  we  see  that  for  a  fixed  number  of  transmit  antennas 
iVt,  even  as  little  as  two  antennas  per-RRU  can  dramatically  improve  feasibility. 

2.3.4  Algorithm  for  IA  with  Strict  Power  Constraints 

While  Proposition  3  gives  properness  conditions  under  which  IA  with  strict  per-RRU  power 
constraints  is  expected  to  be  possible,  it  provides  no  insight  on  how  to  realize  such  alignment 
precoders.  To  better  demonstrate  the  feasibility  of  I A  with  per-RRU  constraints,  we  extend 
the  method  of  alternating  minimization  used  in  [51,64]  to  provide  a  simple  algorithm  which 
satisfies  both  the  alignment  conditions  in  (35)  and  the  per-RRU  constraints  in  (42). 

The  alternating  minimization  strategy  used  in  [51,64]  can  be  summarized  as  iteratively 
minimizing  the  total  power  of  leakage  interference,  defined  as 

Jil  ({FJt€K  ,  {W,}ler)  =  Y,  E  llw^H«F'llt.  (48) 

eeic  keK\£ 

over  alternating  choices  of  W*,  and  F*..  We  refer  the  reader  to  [51,64]  for  a  detailed  derivation 
of  the  optimal  choice  of  W k  and  in  each  iteration  and  give  the  final  result  here  for  brevity. 
At  each  iteration,  the  combiners  W*,  are  chosen  as 
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followed  by  an  update  to  the  choice  of  precoders  given  by 


Ffc  =  uN? 

K  mm 


(50) 


where  ( • )  denotes  the  Ns  least  dominant  eigenvectors  of  a  symmetric  matrix.  To  satisfy 
the  per-RRU  power  constraints,  we  add  a  projection  step  onto  the  set  of  precoders  with 
a  fixed  RRU  transmit  power  of  .  The  projection  on  that  set  amounts  to  scaling  and 
updating  the  blocks  of  as 


ih) 


<(— 


•Wrru  p(r) 


p(r)| 


The  steps  of  the  algorithm  are  listed  more  formally  in  Algorithm  1. 
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2.3.5  Intuition  Obtained  from  the  Results 

Our  results  indicate  that  distributed  antennas  are  an  attractive  system-level  solution  to  in¬ 
crease  the  utility  of  IA  in  practical  wireless  systems.  While  maximum  power  constraints  may 
reduce  sum-rate  relative  to  unconstrained  IA,  and  strict  power  constraints  may  reduce  feasi¬ 
bility,  the  system  configurations  we  considered  showed  promising  sum-rate  performance  for 
most  users  in  a  wireless  network.  As  a  result,  the  SNR  boost  provided  by  DAS  is  sufficient 
to  justify  the  “reduced”  sum-rate  and  feasibility.  Further,  we  showed  that  existing  IA  algo¬ 
rithms  can  be  extended  to  systems  with  per-radio  power  constraints,  thus  the  performance 
demonstrated  in  theory  can  be  achieved  by  practical  iterative  algorithms.  In  conclusion,  our 
work  proves  that  the  combination  of  IA  and  DAS  can  improve  data  rates  for  the  major¬ 
ity  of  mobile  users  when  compared  to  (i)  IA  in  systems  with  co-located  antennas,  and  (ii) 
distributed  antennas  with  existing  interference  management  solutions. 


2.4  Prototyping  MIMO  Interference  Alignment 

Most  prior  work  on  IA  makes  assumptions  about  perfect  synchronization  and  CSI  at  the 
transmitter  (CSIT),  which  is  not  realistic.  To  overcome  these  issues,  we  developed  a  proto¬ 
type  of  interference  alignment  that  operates  in  a  completely  distributed  fashion,  with  each 
node  performing  estimation  and  synchronization. 

CSI  feedback  quality  was  studied  with  limited  feedback  [66]  and  analog  feedback  [37] .  The 
impact  of  feedback  combined  with  imperfect  channel  estimation  was  established  in  [70].  In 
[70],  it  is  assumed  that  the  CSI  feedback  is  given  with  analog  feedback,  and  IA’s  performance 
with  imperfect  feedback  is  analyzed  with  the  CSI  errors  caused  by  the  channel  estimation 
error  in  feedforward  channel,  the  channel  estimation  error  in  feedback  channel,  and  the 
decoding  error  of  CSI  data  in  feedback  channel.  A  main  conclusion  is  that  analog  feedback 
is  a  good  option  for  implementing  interference  alignment.  Consequently  we  incorporate 
analog  interference  alignment  into  our  prototype. 

There  has  been  limited  work  on  IA  from  an  experimental  perspective,  essentially  reported 
in  [71-73]  and  [11],  IA  prototypes  are  built  and  the  performance  of  IA  is  measured  and  studied 
from  them.  In  [71],  the  feasible  setups  for  IA  for  wireless  LAN  system  and  the  imperfect 
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practical  issues  such  as  synchronization  errors  are  studied.  In  [72],  the  wireless  channel  is 
measured  from  the  authors’  prototype,  and  the  sum  rate  is  calculated  from  it  mainly  showing 
that  the  gain  in  sum  rate  achieves  under  the  measured  channel.  In  [73] ,  including  our  initial 
results  reported  in  [If],  the  CSI  is  measured  at  the  receivers,  precoding  vectors  are  obtained 
with  the  CSI,  and  the  sum  rate  is  calculated  from  over-the-air  precoded  transmission. 

The  main  contribution  of  our  work  is  the  first  implementation  and  measurement  of  real¬ 
time  distributed  IA  system.  Though  there  has  been  prior  experimental  work,  none  has 
considered  the  practical  issues  when  the  nodes  of  I A  are  physically  distributed,  and  con¬ 
sequently  the  system  experiences  overhead  and  loss  of  accuracy  when  sharing  CSI.  In  our 
prototype,  we  implemented  over-the-air  time  and  frequency  synchronization  and  feedback 
mechanisms  to  achieve  this  goal.  From  the  prototype,  the  sum  rate  was  measured  showing 
that  IA  has  a  significant  gain  over  the  issues.  Other  measurements  were  performed  con¬ 
sidering  the  accuracy  of  IA  solution,  synchronization  and  CSI  feedback  mainly  showing  the 
relationship  between  overhead  and  performance. 

2.4.1  Mathematical  Framework  and  Problem  Statement 

Our  system  models  the  K-user  interference  channel.  For  the  interference  channel,  there  is  one 
by  one  mapping  of  a  transmitter  and  a  receiver.  Only  the  data  streams  from  i-th  transmitter 
is  the  desired  streams  to  i-th  receiver,  and  the  streams  from  the  other  transmitters  are  the 
interfering  streams.  Each  of  the  transmitters  and  the  receivers  has  Nt  and  Nr  antennas 
respectively,  and  all  nodes  are  separated.  Our  hardware  uses  Nt  =  Nr  =  2. 

IA  basically  assumes  perfect  synchronization  and  perfect  CSI  measurement  and  feedback, 
which  are  not  feasible  in  real  world.  The  requirement  for  time  and  frequency  synchronization 
for  IA  is  studied  in  [71].  It  is  known  from  [71]  that  the  phase  offsets  caused  by  time  and 
frequency  synchronization  errors  are  not  a  problem  for  IA.  The  reason  is  that  the  phase 
error  in  synchronization  does  not  affect  the  antenna  subspaces  where  I A  works.  Besides  this 
requirement  in  synchronization,  a  protocol  for  over-the-air  nodes’  synchronization  is  required 
for  the  distributed  system. 

There  are  two  types  of  synchronization  for  the  distributed  IA  system:  synchronization 
among  transmitters  and  synchronization  between  the  transmitter  side  and  the  receiver  side. 
Once  the  former  is  achieved,  then  the  latter  becomes  a  point  to  point  synchronization  that 
is  well  studied  in  literature.  Thus,  most  of  the  effort  in  our  prototype  was  spent  on  synchro¬ 
nization  among  the  transmitters. 

For  our  interference  channel  model,  all  the  transmitters  use  the  same  time  and  frequency 
resources  to  send  their  own  streams,  so  the  multiple  streams  are  overlapped  at  the  receivers. 
If  the  transmitters  are  asynchronous,  then  the  receivers  have  to  estimate  and  compensate 
for  different  time  and  frequency  offsets  for  each  of  the  transmitted  streams,  and  it  is  an 
expensive  task  for  the  receivers  as  an  example  of  asynchronous  multiuser  MIMO  system  in 
[74],  Otherwise,  if  the  transmitters  are  globally  synchronized,  the  receivers  need  to  estimate 
and  compensate  only  the  global  time  and  frequency  offsets  for  all  the  overlapped  streams, 
and  it  has  lower  synchronization  overhead  than  the  case  of  the  asynchronous  transmitters. 

For  the  transmitter  synchronization,  master-slave  synchronization  protocols  have  been 
proposed  in  [75,76]  and  [77].  In  [75],  a  precise  symbol-level  time  synchronization  protocol  is 
introduced  for  wireless  LAN  systems,  and  in  [76]  a  frequency  synchronization  method  based 
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Dj :  Propagation  delay  between  the  master  TX  and  the  slave  TX  i 
Tp,nm :  Propagation  delay  between  RX  n  and  TX  m 


(a)  System  with  different  propagatin  delays. 


TX  sends  probe 
RX  receives  probe 
RX  detects  probe 
RX  transmits  feedback 
TX  receives  feedback 
TX  detects  feedback 


(b)  Acquisition  of  one  propagation  delay. 


Figure  10:  Over-the-air  master-slave  synchronization  protocol. 


on  the  time  synchronization  protocol  in  [75]  is  also  presented.  Fig.  10(a)  shows  an  example 
of  the  system  that  has  different  propagation  delays  between  any  two  nodes.  The  protocol 
in  [75]  assumes  different  propagation  delays.  One  of  the  transmitters  works  as  the  master 
transmitter  while  others  become  the  slave  transmitters.  The  goal  is  to  align  the  transmitters 
in  time  and  frequency  so  that  the  signals  from  all  the  transmitters  arrive  as  close  as  possible 
in  time  and  frequency  at  the  receivers. 

To  acquire  the  propagation  delays  and  the  other  internal  delays,  each  of  the  transmitters 
periodically  sends  a  probe  to  the  others.  Fig.  10(b)  shows  the  possible  delay  elements  and  the 
method  how  to  obtain  them.  Once  transmitter  sends  the  probe  (this  transmitter  is  denoted 
as  ’TX’  in  Fig.  10(b)),  and  then  each  of  the  other  transmitters  and  the  receivers  (denoted 
as  ’RX’  in  Fig.  10(b))  senses  the  probe,  measures  the  internal  processing  delay  to  detect 
the  probe  (T§)  and  the  turnaround  delay  from  receiving  mode  to  transmitting  mode  (TV), 
informs  the  sum  of  those  delays  back  to  the  transmitter  that  sends  the  probe.  Then,  the 
transmitter  receives  this  feedback  from  each  of  the  nodes,  and  measures  the  air  propagation 
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delays  (Tp)  with  its  internal  detection  time  (Tp)  and  the  delay  information  (T§  +  T't)  given 
by  the  feedback  node.  The  other  transmitters  perform  the  same  task  in  turn. 

When  all  the  required  delay  information  is  obtained,  the  system  is  ready  to  work  in  syn¬ 
chronous  manner.  The  master  transmitter  sends  its  training  to  the  other  transmitters,  and 
each  of  the  other  transmitters  adjusts  the  time  with  the  previously  measured  delay  infor¬ 
mation  to  send  its  own  packet.  After  some  time  margin  to  this  adjustment,  all  transmitters 
send  their  synchronized  packets  to  the  receivers.  If  multiple  receivers  exist  in  the  system, 
however,  it  is  impossible  to  make  all  receivers  receive  the  transmitters’  signals  at  the  exactly 
same  time.  Instead,  the  transmitters  adjust  their  packet  transmission  time  so  that  the  sum 
timing  error  at  all  the  receivers  is  minimized,  e.g.  all  the  signals  are  to  be  received  within 
the  cyclic  prefix  of  orthogonal  frequency  division  multiplexing  (OFDM)  system. 

For  frequency  synchronization,  the  slave  transmitters  measure  the  frequency  offset  with 
the  master’s  training,  recover  the  offset  before  they  send  their  data  packet  [76].  Along  with 
the  frequency  offset,  the  phase  offset  between  the  master  and  the  slaves  are  also  considered 
in  [76].  The  phase  offset  caused  by  oscillator  offset  is  inevitable,  and  it  may  affect  the 
performance  of  wireless  systems.  To  enable  the  slaves  to  measure  the  phase  difference  from 
the  master,  two  identical  trainings  are  sent  from  the  master  :  one  at  the  first  transmission 
of  the  synchronization  training  and  the  other  at  actual  data  transmission.  The  slaves  detect 
these  two  trainings  and  compare  them  to  estimate  phase  offset  between  two.  If  the  residual 
frequency  offset  is  small,  it  can  be  assumed  that  there  is  only  the  constant  phase  offset 
between  the  master  and  a  slave.  The  slaves  also  compensate  the  phase  offset  before  their 
data  transmission. 

In  [75,76],  TDD  is  assumed,  but  there  is  also  a  scheme  that  assumes  frequency  division 
multiplexing  (FDD)  and  continuous  transmission  [77].  In  [77],  the  master  continuously 
sends  its  OFDM-based  training  to  the  slaves  via  a  dedicated  synchronization  channel,  and 
the  phase  rotation  of  OFDM  subcarriers  caused  by  time  and  frequency  offset  from  the  master 
is  measured  and  compensated  by  the  slaves  before  they  send  their  own  data  stream. 

For  CSI  feedback,  there  are  three  major  error  types:  the  error  caused  by  AWGN,  the 
limited  feedback  error,  and  the  delayed  feedback  error.  The  limited  feedback  error  occurs 
because  of  the  quantization  of  the  estimated  CSI,  while  the  delayed  feedback  error  is  by 
channel’s  time- varying  property.  IA  with  limited  feedback  is  first  analyzed  in  [66].  In  [66], 
the  CSI  quantization  error  with  Grassmannian  codebook  [78]  is  analyzed  assuming  perfect 
CSI  measurement  and  feedback.  According  to  its  result,  the  DOF  of  K/2  for  single  input 
and  single  output  IA  system  is  achievable  if  the  number  of  total  feedback  bits  is  larger  than 
K(L  —  l)logPf,  where  L  is  the  number  of  multipaths  in  feedforward  channel  and  Pf  is  the 
trasnmit  power  of  CSI  feedback. 

Our  prior  work  on  analog  feedback  with  IA  is  presented  in  [37]  as  a  way  to  flexibly 
implement  feedback  in  a  way  that  avoids  quantization  error.  In  [37],  the  feedback  channel’s 
estimation  error  by  AWGN  at  feedback  training  and  the  feedback  data  error  caused  by 
AWGN  at  data  symbols  are  considered.  It  is  shown  that  the  multiplexing  gain  of  IA  is 
preserved  only  with  a  constant  sum  rate  loss  when  the  transmitting  power  of  feedback  linearly 
relates  to  the  transmitting  power  of  feedforward.  The  upper  bound  of  the  constant  sum 
rate  loss  A Rsum  is  given  as  the  function  of  the  length  of  the  training  for  feedback  channel 
estimation  (rp),  the  length  of  the  CSI  feedback  (rc),  and  the  ratio  of  Pf  and  the  transmit 
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Figure  11:  Operation  scenario  for  the  prototype. 


power  P  of  feedfoward  channel: 

_  p  2 

yy/o^l  +  —c(tp,tc)(K  -  — )).  (52) 

-*  /  di 

%  j 

where  c(rp,  rc)  is  a  value  that  depends  on  rp  and  rc,  and  if  either  rp  or  rc  gets  larger,  c(tp,  tc) 
gets  smaller  showing  that  with  more  overhead  in  the  feedback  channel  estimation  training 
and  feedback  data,  the  sum  rate  loss  A Rsum  gets  smaller.  To  maintain  the  multiplexing 
gain,  Pf  needs  to  linearly  scale  with  P ,  i.e.  Pf  =  aP.  Otherwise,  if  Pf  is  not  linear  with 
P  as  Pf  =  aPf  there  exists  a  loss  in  multiplexing  gain  according  to  ft.  A  similar  extended 
study  also  considering  the  feedforward  channel  estimation  error  is  presented  in  [70]. 

2.4.2  System  Realization 

We  are  targeting  2x2  MIMO  three  user  IA  system.  Nt  =  Nr  =  2  for  all  user  pairs,  Ns  =  1, 
and  K  =  3.  OFDM  is  used  for  both  feedforward  and  feedback  channels.  With  this  setup,  two 
antenna  subspaces  are  available  for  a  user  pair,  and  two  interfering  streams  are  aligned  in  one 
subspace  at  each  receiver  while  the  desired  stream  is  received  in  the  other  subspace.  Both 
the  closed  form  solution  and  the  iterative  method  in  [79]  are  implemented  in  our  prototype. 

There  are  three  main  phases  in  IA  operation  scenario:  training,  feedback  and  data  phases. 
Fig.  11  summarizes  the  signals  that  are  needed  to  be  transmitted  at  each  of  the  phases,  and 
also  shows  the  tasks  that  should  be  done  during  each  of  the  phases  at  the  transmitters  and 
the  receivers. 

•  Training  Phase:  This  is  mainly  for  CSI  measurement  at  the  receivers.  The  syn¬ 
chronization  training  is  sent  from  the  transmitters  followed  by  the  CSI  training.  The 
receivers  first  synchronize  to  the  transmitters  in  time  and  frequency  with  the  received 
synchronization  training,  and  then  measure  CSI  which  should  be  given  back  to  the 
transmitters  at  the  feedback  phase.  Furthermore,  the  synchronization  among  the  trans¬ 
mitters  should  be  done  at  this  phase.  The  details  about  the  training  structure  and  the 
transmitter  synchronization  is  given  in  the  following  section. 

•  Feedback  Phase:  This  phase  is  for  CSI  feedback.  Each  of  the  receivers  makes  a 
packet  which  is  consisted  of  the  synchronization  training,  the  CSI  training,  and  the  CSI 
feedback  data,  and  the  packet  is  sent  to  the  transmitters  via  over-the-air  transmission. 
The  packets  are  sent  in  time  orthogonal  manner.  After  receiving  and  decoding  these 
packets,  the  transmitters  calculate  IA  precoding  vectors. 
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Figure  12:  Transmitter  synchronization  in  our  prototype. 

•  Data  Phase:  This  phase  is  for  the  precoded  training  and  data  transmission.  The 
synchronization  training,  the  precoded  CSI  training,  and  the  precoded  data  are  sent 
to  the  receivers  in  order.  The  receivers  first  synchronize  to  the  transmitters,  calculate 
the  decoding  vectors  with  the  precoded  CSI  training,  and  decode  the  precoded  data 
with  them.  Finally,  the  sum  rate  is  calculated  and  the  operation  ends. 

The  master-slave  protocol  is  used  for  over-the-air  transmitter  synchronization  in  our 
prototype,  but  there  are  two  differences  from  the  algorithms  in  [75]  and  [76].  First  of  all, 
the  phase  synchronization  in  [76]  is  not  required  in  our  prototype  since  it  does  not  affect  IA 
performance.  Also,  time  synchronization  is  simplified  from  [75].  The  over-the-air  propagation 
delays  between  nodes  are  assumed  to  be  the  same.  With  this  assumption,  there  is  no  need  to 
measure  the  air  propagation  delays.  Then,  the  processing  delay  and  the  turn-around  delay  at 
both  the  transmitters  and  the  receivers  do  not  matter  only  if  the  following  two  requirements 
are  guaranteed  :  1)  a  symbol  should  be  sent  from  transmitting  antennas  exactly  at  desired 
time  at  the  transmitters,  and  2)  the  exact  time  when  a  symbol  arrives  to  receiving  antennas  is 
known  at  the  receivers.  These  requirements  are  supported  by  our  implementation  hardware 
and  software. 

Fig.  12  illustrates  our  protocol.  For  our  distributed  system,  it  is  assumed  that  each  of 
the  nodes  are  physically  separated,  and  this  means  that  every  node  works  with  its  own  time 
and  frequency  references  not  knowing  the  others’  references.  At  the  training  phase  which  is 
the  starting  point  of  one  turn  of  IA  procedure,  the  synchronization  is  done  in  the  following 
three  steps: 

•  Step  1.  The  master  transmitter  first  sends  its  training  to  the  slave  transmitters  at  time 
zero  with  its  local  timer.  When  each  of  the  slaves  receives  and  detects  the  training, 
it  knows  its  local  time  Tj^  when  the  first  sample  of  the  training  is  arrived  at  its 
antennas.  Here,  n  is  the  slave  index.  The  slaves  also  measure  the  frequency  offset  from 
the  master’s  frequency  with  this  training. 

•  Step  2.  Tq  is  the  waiting  time  for  synchronized  transmission,  and  it  is  known  to  all 
transmitters.  The  master  and  the  slaves  wait  until  their  local  time  reaches  Tc  for  the 
master,  and  +  Tc  for  the  slaves. 
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•  Step  3.  All  the  transmitters  send  their  own  training  to  the  receivers  immediately  when 
their  waiting  time  ends.  Each  of  the  slaves  recovers  the  measured  frequency  offset  from 
its  own  training  before  it  is  sent.  The  training  and  data  are  now  synchronized  in  time 
and  frequency. 

The  same  operation  is  performed  at  the  data  phase  since  the  time  between  the  training 
phase  and  the  data  phase  may  not  short  enough  for  the  slaves  to  maintain  their  previous 
synchronization.  At  the  data  phase,  the  starting  time  of  the  data  phase,  TD,  is  only  known 
to  the  master  and  it  sends  its  training  again  to  the  slaves  at  that  time.  The  slaves  now  detect 
it  at  time  ,  and  all  transmitters  send  the  synchronization  training  and  data  after 

Tc  time.  Besides  the  transmitter  synchronization,  the  receivers  also  need  to  synchronize  to 
the  transmitters.  Each  receiver  synchronizes  independently  to  the  transmitters,  and  there  is 
no  cooperative  synchronization  among  the  receivers.  The  auto-correlation  method  in  [80]  is 
used  for  synchronization  at  all  nodes. 

The  main  feedback  scheme  in  this  paper  is  analog  feedback.  Under  the  assumption  of 
flat  fading  channel,  receiver  i  measures 

H ij,  j  =  1,2,3.  (53) 

Since  each  H ^  is  2  x  2  matrix,  there  are  K  x  Nt  x  Nr  complex  values  for  the  measured  CSI  at 
each  receiver.  These  values  are  mapped  to  Nsc  subcarriers  by  Nsc  x  (. K  x  Nt  x  Nr)  mapping 
matrix,  and  the  Nsc  subcarriers  form  an  OFDM  symbol.  As  a  result,  each  receiver  has  an 
OFDM  symbol  for  analog  feedback.  The  transmitters  that  receive  the  OFDM  symbols  apply 
the  pseudo  inverse  of  the  mapping  matrix  to  the  received  OFDM  symbols,  and  find  the  CSI 
values.  The  transmit  power  which  decides  the  difference  in  SNR  between  feedforward  and 
feedback  channel  as  in  (52)  is  a  design  parameter. 

For  limited  feedback  as  a  competitor,  the  explicit  feedback  beamforming  method  from 
the  beamforming  specification  of  802. lln  wireless  LAN  is  used  [81].  With  the  quantization 
method,  Ngain  bits  are  used  for  the  amplitude  of  CSI,  and  Nq  bits  are  used  to  quantize  each 
of  real  and  imaginary  values  of  the  measured  CSI.  Ngain  =  3  and  Nq  =  4,  5,  6  and  8  for  our 
implementation.  Gref  is  chosen  to  minimize  the  mean  squared  error  (MSE)  by  quantization. 
Note  that  this  is  really  a  brute  force  quantization  technique  and  is  not  “limited”  in  any  sense. 

2.4.3  Experiments  and  Results 


Table  1:  OFDM  parameters  for  the  prototype. 


FFT  length 

128 

Cyclic  Prefix  Length 

32 

Number  of  null  subcarriers 

23 

Number  of  data  subcarriers 

105 

The  system  is  designed  to  operate  at  2.4  GHz,  the  bandwidth  would  ideally  be  large  as 
20  MHz,  but  it  is  lower  as  1  MHz  in  our  prototype  because  of  hardware  constraints.  Both 
feedforward  and  feedback  channel  use  the  same  frequency,  and  the  system  is  a  TDD  system. 
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Figure  13:  Packet  structures. 


OFDM  is  used  for  the  training  and  data,  and  only  the  synchronization  training  is  generated 
at  time  domain.  The  OFDM  parameters  are  summarized  in  Table  1. 

As  described  in  the  previous  sections,  our  system  works  in  three  phases  :  the  training, 
the  feedback  and  the  data  phases.  For  each  of  the  phases,  the  packet  structure  is  shown  in 
Fig.  13.  In  Fig.  13,  the  packet  structures  for  the  training  phase  (Fig.  13(a)),  the  feedback 
phase  (Fig.  13(b)),  and  the  data  phase  (Fig.  13(c))  are  given.  There  are  three  types  of 
trainings  in  our  system: 

•  Synchronization  training:  The  synchronization  training  is  for  the  nodes’  synchro¬ 
nization.  It  has  two  parts  :  short  training  for  coarse  time  synchronization  and  long 
training  for  fine  time  and  frequency  synchronization.  A  length  17  Zadoff-Chu  sequence 
is  repeated  five  times  as  the  short  training,  and  a  length  29  Zaduff-Chu  sequence  is 
repeated  three  times  as  the  long  training.  The  training  has  repeated  patterns  in  time 
domain,  and  the  synchronization  is  done  by  the  auto-correlation  method  in  [80]. 

•  CSI  training:  The  CSI  training  is  for  CSI  measurement  to  be  sent  back  to  the  trans¬ 
mitters  for  IA  precoding  vector  calculation.  One  CSI  training  is  an  OFDM  symbol, 
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Figure  14:  Software  and  hardware  configuration  of  the  prototype. 

and  it  is  sent  from  each  of  the  transmitting  antennas  in  time  orthogonal  manner  as 
in  Fig.  13(a).  As  such,  six  OFDM  symbols  are  transmitted  for  our  three  user  2x2 
MIMO  system.  The  training  does  not  experience  precoding  nor  decoding.  Thus,  the 
receivers  can  find  the  pure  condition  of  current  wireless  channel.  With  this  training, 
each  receiver  measures  2x6  CSI  matrix. 

•  Precoded  training:  The  precoded  training  is  sent  at  the  data  phase.  It  experiences 
I A  precoding  at  the  transmitters,  where  each  user  in  our  three  user  system  has  one 
data  stream  that  is  mapped  to  two  antennas  via  the  precoding  vector.  It  is  used  for 
the  decoding  vector  calculation.  Also,  since  this  training  is  precoded  in  the  same  way 
with  the  precoded  data,  it  is  used  to  equalize  the  data  symbols.  The  training  is  also 
an  OFDM  symbol,  and  it  is  sent  from  two  transmitting  antennas  of  a  transmitter  at  a 
time.  Each  transmitters  sends  it  in  time  orthogonal  manner  as  in  Fig.  13(c). 

Fig.  14  shows  the  hardware  and  software  configuration  of  our  prototype.  The  system  uses 
two  computers  that  are  equipped  with  dual-core  Intel  Xeon  2.67  GHz  processors  and  12GB  of 
memory.  One  computer  is  used  to  control  the  transmitters,  and  the  other  is  for  the  receivers. 
National  Instruments  (NI)  USRP-2921  [82]  is  used  for  analog  to  digital  convertor  (ADC)  / 
digital  to  analog  convertor  (DAC)  and  radio  frequency  (RF)  front-end.  One  USRP  works 
as  one  antenna,  so  two  USRPs  are  used  for  a  node,  and  accordingly  one  computer  controls 
six  USRPs.  A  USRP  can  switch  between  transmitting  mode  and  receiving  mode,  e.g.  the 
USRPs  for  the  transmitters  work  in  transmitting  mode  for  the  training  and  data  phases,  but 
in  receiving  mode  for  the  feedback  phase.  The  computer  and  the  USRPs  are  connected  via 
TCP/IP  connections.  National  Instruments  (NI)  Lab  VIEW  which  is  a  model-based  digital 
signal  processing  implementation  software  [83]  is  used  for  our  software  defined  radio  (SDR) 
implementation.  The  control  of  USRPs  is  also  done  by  Lab  VIEW’S  USRP  driver.  Though 
only  one  computer  is  used  for  all  the  transmitters,  and  similarly  only  one  computer  is  used 
for  all  the  receivers,  the  signal  processing  for  each  node  is  completely  separated  from  each 
other  not  to  violate  the  interference  channel  model,  i.e.  there  are  three  parallel  independent 
processing  chains  in  Lab  VIEW  on  a  computer. 

Since  our  system  is  distributed  system,  each  node  should  work  with  its  own  oscillator. 


Figure  15:  Hardware  configuration  of  one  node. 


V 


Figure  16:  Test  environment  of  distributed  nodes. 


The  GPS  disciplined  oscillator  (GPSDO)  module  that  can  be  attached  to  USRP  is  used  as 
the  oscillator  for  a  node  [84] .  The  GPSDO  is  originally  for  the  global  synchronization  to  the 
GPS  signal  which  is  only  available  at  outdoors.  There  is,  however,  no  assumption  of  outdoor 
operation  for  IA  system  in  general,  so  the  GPSDOs  are  not  used  to  synchronize  to  the  global 
GPS  signal.  Instead,  a  GPSDO  is  used  to  provide  time  and  frequency  references  only  to  a 
single  node  in  our  prototype.  Though  the  GPSDOs  are  used,  our  distributed  synchroniza¬ 
tion  algorithm  is  designed  assuming  that  each  of  the  GPSDOs  generates  the  independent 
time  reference  (pulse  per  second  (PPS))  and  frequency  reference  (10MHz)  without  being 
synchronized  to  the  GPS  signal  since  our  prototype  is  tested  in  indoor  environment  where 
the  global  GPS  signal  is  unavailable.  The  PPS  and  10MHz  are  the  required  time  and  fre¬ 
quency  references  for  the  USRPs  respectively.  Six  GPSDOs  are  used  because  there  are  six 
nodes  in  the  system,  and  the  time  and  frequency  reference  signals  from  a  GPSDO  are  only 
shared  between  two  USRPs  which  form  a  node’s  two  antennas.  Fig.  15  shows  the  hardware 
configuration  of  a  node  in  the  prototype.  The  MIMO  cable  between  the  USRPs  in  Fig.  15 
is  only  used  for  the  ethernet  connection  of  two  USRPs  to  a  computer. 

Fig.  16  shows  the  schematic  of  distributed  test  environment.  Both  the  transmitters  and 
the  receivers  are  physically  separated  to  verify  our  system.  To  make  the  SNRs  of  the  air 
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(b)  Comparison  of  measured  sum  rate  to  simulated  sum 
rate. 


Figure  17:  IA  performance  with  analog  feedback 


links  be  similar  as  possible,  the  distances  between  any  transmitter  and  any  receiver  are  made 
to  be  the  same  as  possible.  Also,  the  AWGN  generator  which  is  used  to  scale  the  overall 
SNR  of  the  system  is  located  at  the  center  of  the  nodes. 

In  this  paper,  the  sum  rate  of  distributed  IA  with  analog  feedback  is  first  presented. 
Then,  more  detailed  measurements  on  the  relationship  between  the  accuracy  of  IA  solution, 
residual  frequency  offset,  and  CSI  feedback  are  also  given.  For  all  the  results  in  this  paper, 
the  distributed  nodes’  synchronization  protocol  is  applied. 

The  sum  rate  of  our  system  is  given  in  Fig.  17.  Each  dot  in  the  first  subfigure  is  the 
measured  instantaneous  sum  rate  within  SNR  range  from  15dB  to  35dB  and  the  solid  line 
shows  their  linear  fitting,  and  the  second  subfigure  shows  the  comparison  of  the  linear  fitting 
to  the  simulation  results.  The  simulation  is  performed  under  different  setups  :  1)  perfect 
CSI  estimation  (CE)  and  feedback,  2)  only  with  feedforward  CE  error,  3)  with  feedforward 
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Figure  18:  Sum  rate  with  iterative  IA  method. 

CE  and  feedback  data  errors,  4)  with  the  errors  in  3)  and  feedback  CE  error,  and  5)  with 
all  the  error  sources  :  feedforward  CE  error,  feedback  data  error,  and  feedback  CE  error, 
and  feedforward  precoded  CE  error.  It  is  observed  from  Fig.  17  that  the  multiplexing  gain 
is  maintained  by  analog  feedback  only  with  a  constant  sum  rate  degradation  if  the  feedback 
transmit  power  is  linear  to  the  feedforward  transmit  power  for  CSI  measurement  as  analyzed 
in  [37].  This  result  is  possible  since  the  SNRs  for  feedforward  channel  and  feedback  channel 
are  the  same  in  our  experiments. 

To  study  the  performance  with  imperfect  IA  solution,  Fig.  18  shows  the  sum  rate  versus 
the  number  of  iterations  for  the  iterative  IA  method.  The  iterative  method  from  [79]  is 
used,  and  the  number  of  iterations  is  an  important  design  parameter  for  the  method.  The 
mean  sum  rate  at  SNR  =  30dB  and  SNR  =  20dB  are  measured  and  compared  to  those  from 
the  simulation  with  the  perfect  CSI  feedback.  Though  the  sum  rate  from  the  prototype  is 
smaller  than  that  from  the  simulation  because  it  is  limited  by  CSI  estimation  and  feedback, 
the  number  of  iterations  to  achieve  the  sum  rate  is  not  reduced.  This  shows  that  the  leakage 
from  the  imperfect  IA  solution  is  independent  to  the  other  error  sources  regardless  of  SNR 
of  the  system. 

Fig.  19  shows  the  IA  performance  with  synchronization  error.  The  time  synchronization 
error  of  IA  with  OFDM  does  not  affect  the  IA  performance  only  if  the  timing  synchroniza¬ 
tion  is  done  only  within  cyclic  prefix,  and  this  requirement  is  not  new  to  IA.  Frequency 
synchronization,  however,  can  not  be  perfect  in  the  real  systems,  and  there  always  remains 
the  residual  frequency  offset  after  synchronization.  It  causes  inter  channel  interference  (ICI) 
which  can  not  be  recovered  by  IA,  so  it  is  added  as  an  extra  error  source  to  the  system. 
Furthermore,  the  ICI  by  the  residual  frequency  offset  affects  all  the  CSI  measurement  and 
feedback  steps  which  are  listed  in  Fig.  17.  In  Fig.  19,  the  sum  rate  is  plotted  with  different 
%  frequency  offsets  with  respect  to  the  subcarrier  spacing.  The  same  frequency  offset  is 
added  to  the  air  links  between  the  master  transmitter  and  the  slave  transmitters,  and  to  the 
connection  between  the  transmitters  and  the  receivers.  The  SNR  is  30dB. 

The  analog  feedback  is  also  compared  to  limited  feedback.  For  limited  feedback,  the 
uncompressed  explicit  scalar  quantization  in  802. lln  specification  [81]  is  used  as  described 
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Figure  19:  IA  performance  with  the  residual  frequency  synchronization  error 


in  the  previous  section.  Nq  =  4,5,6  and  8  are  applied  for  the  evaluation.  It  is  observed  from 
Fig.  20(a)  that  the  slope  of  sum  rate  increases  as  more  bits  are  allocated,  i.e.  as  feedback 
overhead  is  increased,  which  means  more  bits  are  needed  as  feedback  as  SNR  gets  high  to  get 
the  multiplexing  gain  of  IA.  Fig.  20(b)  shows  the  sum  rate  against  MSE  of  limited  feedback. 
The  fitting  line  is  the  linear  connection  of  the  means  of  the  sum  rates  with  different  Nq.  The 
results  in  Fig.  20  gives  us  some  interesting  observations.  First  of  all,  the  multiplexing  gain 
of  quantization-based  feedback  increases  with  increasing  number  of  Nq,  which  corresponds 
to  the  previous  analysis  in  [66] .  Furthermore,  compared  to  the  quantization-based  feedback, 
analog  feedback  has  higher  sum  rate  performance  at  high  SNR  region  than  the  quantization- 
based  feedback.  This  is  because  the  MSE  of  analog  feedback  is  reduced  with  increasing 
feedback  SNR,  but  the  MSE  of  quantization-based  feedback  is  fixed  regardless  of  SNR  if  it 
is  assumed  that  there  is  no  bit  error  for  quantization-based  feedback. 

2.4.4  Intuition  Obtained  from  the  Results 

We  implemented  a  prototype  of  a  three  user  IA  system  with  2x2  antenna  configuration  for 
each  user.  The  prototype  is  fully  distributed  and  uses  wireless  both  for  direct  communication 
and  feedback.  With  the  prototype,  the  performance  of  IA  with  the  practical  issues  was 
studied,  with  an  emphasis  on  the  feedback  quality.  The  nodes  were  assumed  to  be  physically 
separated  and  to  work  independently  without  cooperation  following  the  basic  assumption  of 
interference  channel.  The  distributed  operations  were  achieved  by  developing  over-the-air 
time  and  frequency  synchronization  protocol  and  feedback  methods.  Our  results  show  that 
it  is  possible  to  implement  IA  in  a  distributed  fashion.  Further,  they  show  the  efficacy  of 
analog  feedback  in  a  real-world  setting.  This  means  it  is  a  viable  approach  for  achieving  good 
performance  as  an  alternative  to  quantization  based  approaches.  Finally,  we  confirmed  that 
there  is  a  lot  of  overhead  associated  with  the  protocol  used  for  distributed  synchronization. 
More  work  is  needed  in  the  future  to  design  a  MAC  protocol  that  allows  fast  setup  of 
distributed  IA  clusters. 
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(a)  Measured  sum  rate  with  limited  feedback  and  analog 
feedback 


o 
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(b)  Sum  rate  vs.  MSE  of  limited  feedback. 

Figure  20:  I A  performance  with  limited  feedback 
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